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The  focus  of  this  research  has  been  on  the  assessment  of  the  effects  of  internal  friction  on  the  response  of  aeroelastic  systems  exhibiting 
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major  efforts,  two  computational/theoretical  and  one  experimental: 

(1)  the  analysis  of  the  effects  of  friction  on  a  structural  dynamic  systems  in  which  the  effects  of  the  aerodynamic  forces  have  been  modeled. 
The  computational  expediency  of  this  simple  model  allowed  to  analyze  and  study  a  vast  array  of  cases  from  which  broad  conclusions  were 
derived  the  validity  of  which  was  later  extended  to  the  complex  models  of  task  (2). 

(2)  the  analysis  of  the  effects  of  friction  on  actual  aeroelastic  systems,  i.e.  airfoil  and  flat  plate  in  a  uniform  flow,  in  which  the  aerodynamic 
forces  are  computed  in  time  in  parallel  to  the  structural  dynamic  analysis  in  a  tight  coupling  format. 
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similar  to  the  ones  studied  in  tasks  (1)  and  (2).  This  effort  was  followed  by  a  data  analysis  that  revealed  mostly  similarities,  albeit  a  few 
differences,  with  the  theoretical  results  of  tasks  (1)  and  (2)  above. 
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FINAL  REPORT 


EXECUTIVE  SUMMARY 

Focus  of  the  Present  Investigation 

The  focus  of  this  research  has  been  on  the  assessment  of  the  effects  of  internal 
friction  on  the  response  of  aeroelastic  systems  exhibiting  either  explosive  flutter  or  limit 
cycle  oscillations  to  provide  a  confirmation  of  the  potential  of  friction  as  a  significant 
stabilizing  factor  in  the  limit  cycle  oscillations  observed  on  several  aircraft,  e.g.  F-16. 
The  work  performed  in  this  regard  can  broadly  be  divided  into  three  major  efforts,  two 
computational/theoretical  and  one  experimental: 

(1)  the  analysis  of  the  effects  of  friction  on  a  structural  dynamic  systems  in  which 
the  effects  of  the  aerodynamic  forces  have  been  modeled.  Specifically,  a  dashpot  of 
negative  constant  was  introduced  to  model  the  unstable  linear  aerodynamic  effects,  while 
an  additional  van  der  Pol  restoring  force  was  included  to  characterize  possible  stable 
nonlinear  aerodynamic  effects.  Although  approximate  (because  of  the  modeling  of  the 
aeroelastic  forces),  this  system  was  shown  to  behave  almost  exactly  as  the  aeroelastic 
systems  of  task  (2).  The  computational  expediency  of  this  simple  model  allowed  to 
analyze  and  study  a  vast  array  of  cases  from  which  broad  conclusions  were  derived  the 
validity  of  which  was  later  extended  to  the  complex  models  of  task  (2). 

(2)  the  analysis  of  the  effects  of  friction  on  actual  aeroelastic  systems,  i.e.  airfoil 
and  flat  plate  in  a  uniform  flow,  in  which  the  aerodynamic  forces  are  computed  in  time  in 
parallel  to  the  structural  dynamic  analysis.  A  tight  coupling  between  structural  motions 
and  flow  field  was  achieved  by  proceeding  iteratively  at  each  time  step  until  both  sets  of 
field  equations  (structural  and  aerodynamic)  were  satisfied. 
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(3)  the  experimental  testing  in  the  DLR-  Gottingen  (Germany)  Transonic  wind 
tunnel  of  a  NLR7301  airfoil  equipped  with  a  friction  device  similar  to  the  ones  studied  in 
tasks  (1)  and  (2).  This  effort  was  followed  by  a  data  analysis  that  revealed  mostly 
similarities,  albeit  a  few  differences,  with  the  theoretical  results  of  tasks  (1)  and  (2) 
above. 

The  overall  results  of  this  study  demonstrate  that  friction  can  indeed  provide  a 
stabilization  of  an  impending  flutter  and  can  significantly  decrease  the  amplitude  of 
existing  limit  cycle  oscillations  of  aeroelastic  systems  with  an  appropriate  selection  of  the 
friction  device  parameters  most  notably  natural  frequency  and  coefficient  of  friction. 

Participants  and  Publications 

These  efforts  were  carried  out  by  the  principal  investigators.  Prof.  Danny  D.  Liu 
and  Marc  P.  Mignolet  from  Arizona  State  University  and  the  former  graduate  students 
Anthony  M.  Agelastos  and  Goang  Gae  Choi  both  of  whom  have  successfully  defended 
their  M.S.  thesis  on  this  topic.  The  executive  summary  of  results  presented  below 
highlights  the  findings  discussed  in  details  in  their  theses: 

1.  Agelastos,  A.M.,  Effects  of  Coulomb  Friction  on  Flutter  and  Limit  Cycle  Oscillations 
According  to  a  Structural  Dynamic  Model,  M.S.,  Arizona  State  University,  May  2005. 

2.  Choi,  G.G.,  Effects  of  Coulomb  Friction  on  Aeroelastic  Systems,  M.S.,  Arizona  State 
University,  August  2005. 

and  which  are  included  to  this  final  report  as  Attachments.  Several  additional  persons 
contributed  very  significantly  to  the  experimental  component  of  this  project,  most 
specifically  Dr.  Dallas  Kingsbury  from  Arizona  State  University  and  the  staff  of  the 


DLR-Gottingen  Institute  for  Aeroelasticity,  in  particular  Dr.  Guido  Dietz  and  Dr.  Gunter 
Schewe. 

In  addition  to  the  above  theses,  several  conference  papers  have  been  published 
that  describe  the  findings  of  this  project: 

1.  Choi,  G.G.,  Agelastos,  A.M,  Mignolet,  M.P.,  and  Liu,  D.D.,  “On  the  Impact  of 
Internal  Friction  on  Flutter  Onset  and  Limit  Cycle  Oscillations  Amplitude,”  International 
Forum  on  Aeroelasticity  and  Structural  Dynamics  2005,  Munich,  Germany,  Jun.  28-Jul. 
1,2005.  Paper  IF-022. 

2.  Kingsbury,  D.W.,  Agelastos,  A.M.,  Dietz,  G.,  Mignolet,  M.P.,  Liu,  D.D.,  and  Schewe, 
G.,  “Limit  Cycle  Oscillations  of  Aeroelastic  Systems  with  Internal  Friction  in  the 
Transonic  Domain  -  Experimental  Results,”  Proceedings  of  the  46th  Structures, 
Structural  Dynamics,  and  Materials  Conference,  Austin,  Texas,  Apr.  18-21,  2005. 
AIAA  Paper  AIAA-2005-1914. 

3.  Choi,  G.G.,  Agelastos,  A.M.,  Mignolet,  M.P.,  and  Liu,  D.D.,  “Effects  of  Internal 
Friction  on  the  Dynamic  Behavior  of  Aeroelastic  Systems,”  Proceedings  of  the  45th 
Structures,  Structural  Dynamics,  and  Materials  Conference,  Palm  Springs,  California, 
Apr.  19-22,  2004.  Paper  AIAA-2004- 1591. 

4.  Mignolet,  M.P.,  Agelastos,  A.M.,  and  Liu,  D.D.,  “Impact  of  Frictional  Structural 
Nonlinearity  in  the  Presence  of  Negative  Aerodynamic  Damping,”  Proceedings  of  the 
44th  Structures,  Structural  Dynamics,  and  Materials  Conference,  Norfolk,  Virginia,  Apr. 
7-10, 2003.  Paper  AIAA-2003-1428. 

Finally,  the  following  three  journal  papers  are  in  preparation: 


1.  Kingsbury,  D.W.,  Agelastos,  A.M.,  Dietz,  G.,  Mignolet,  M.P.,  Liu,  D.D.,  and  Schewe, 
G.,  “Measurements  of  Limit  Cycle  Oscillations  of  Aeroelastic  Systems  Induced  by 
Internal  Friction  in  the  Transonic  Domain  -  Experimental  Results”. 

2.  Agelastos,  A.M.,  Choi,  G.G.,  Mignolet,  M.P.,  Liu,  D.D.,  and  Liao,  Y.,  “Effects  of 
Internal  Friction  on  the  Behavior  of  Aeroelastic  Systems.  Part  I:  Increased  Stability  and 
Limit  Cycle  Oscillations  due  to  Friction  only”. 

3.  Choi,  G.G.,  Agelastos,  A.M.,  Mignolet,  M.P.,  and  Liu,  D.D.,  “Effects  of  Internal 
Friction  on  the  Behavior  of  Aeroelastic  Systems.  Part  II:  Reduction  of  Existing  Limit 
Cycle  Oscillations”. 

Summary  of  Findings 

The  findings  of  the  three  separate  tasks  (l)-(3)  will  be  summarized  in  order.  The 
structural  dynamic  model  of  task  (1)  was  a  two  degree  of  freedom  system  in  which  the 
largest  mass  represented  the  wing  while  the  smallest  one  would  model  the  sliding 
component,  e.g.  a  missile  or  store.  The  negative  dashpot  and  the  van  der  Pol  force  were 
assumed  to  act  only  on  the  largest  mass  (wing).  Further,  the  sliding  mass  was  assumed  to 
be  connected  to  the  wing  through  a  friction  slider,- a  spring  ( k j  >0),  and  possibly  a 
dashpot  (positive  viscous  damping,  c\  >0). 

The  system  was  analyzed  first  without  a  van  der  Pol  term.  Further,  in  the  absence 
of  elastic  deformation  during  sliding,  i.e.,  for  Aq  =  0 ,  it  was  shown  that  an  exact 
procedure  can  be  followed  to  transform  the  nonlinear  equations  of  motion  into  a  set  of 
nonlinear  algebraic  equations  and  that  a  stabilization  of  the  unstable  aerodynamic  forces 
was  possible.  The  most  negative  damping  that  could  be  stabilized  in  this  configuration 

was  in  fact  shown  to  be  2rht  (for  r  « 1)  where  r  denotes  the  ratio  of  the  two  masses 


of  the  system.  The  system  was  found  to  exhibit  periodic  motions  with  possible  stick 
phases.  The  introduction  of  a  nonzero  spring  constant  was  shown  to  lead  to  an  increased 
stabilization  potential  that  is  maximum  in  a  “tuned  damper”  configuration,  i.e.  when  the 
natural  frequency  of  the  primary  component  (the  wing  alone)  closely  matches  the  natural 
frequency  of  the  secondary  alone  (sliding  mass  on  spring  /q  ).  In  addition  to  this  benefit, 
a  broader  set  of  motions  was  also  observed.  Besides  periodic  (referred  to  as  single 
frequency  here)  solutions,  multiple  frequency  motions  were  also  found  to  exist, 
especially  near  the  tuned  damper  configuration  while  the  single  frequency  solutions  were 
typically  obtained  near  the  stability  border.  Both  continuous  slip  and  stick  slip  multiple 
frequency  motions  were  noted  which  were  shown  to  be  either  aperiodic  or  chaotic. 
Further,  the  transition  from  single  to  multiple  frequency  was  shown  to  be  a  Hopf 
bifurcation.  It  was  finally  observed  that  the  maximum  steady  state  response  of  the  wing  is 
typically  minimum  for  the  single  frequency  motions  occurring  just  next  to  the 
bifurcation.  A  sharp  increase  in  response  level  is  obtained  after  bifurcation  that  results 
from  the  beating  induced  by  the  presence  of  multiple  frequencies. 

The  inclusion  of  a  van  der  Pol  term  did  not  significantly  change  the  above 
findings  with  the  noted  exception  that  the  system  is  always  stable  even  in  the  absence  of 
friction  in  which  case  it  experiences  a  continuous  slip  single  frequency  motion.  In  fact,  it 
was  observed  that  the  van  der  Pol  term  leads  to  a  significant  reduction  of  the  large 
amplitude  beating  excursions  obtained  otherwise  in  connection  with  multiple  frequency 
solutions  and  appears  to  postpone  the  Hopf  bifurcation.  In  regards  to  stick  slip  vs. 
continuous  slip,  the  decrease  in  response  implied  by  the  presence  of  the  van  der  Pol  term 
was  found  to  increase  the  likelihood  of  stick  slip  motions  which  typically  occur  at  lower 


response  levels  than  continuous  slip  solutions  do.  The  effects  of  the  coefficients  of 
friction  (the  static  and  dynamic  coefficients  were  assumed  equal  here)  was  studied  next. 
It  was  first  noted  that  the  amplitude  of  continuous  slip  motions  (single  or  multiple 
frequency)  is  linearly  scaled  by  the  coefficient  of  friction  and  this  parameter  plays  no 
further  role  in  these  cases.  When  stick  phases  occur,  a  more  complex  dependency  on  the 
coefficient  of  friction  is  obtained.  Next,  it  was  found  that  2  stable  solutions  can  be 
obtained  for  low  coefficients  of  friction.  One  such  solution  has  a  low  amplitude  and  is 
dominated  by  friction  effects  alone  with  little  influence  of  the  van  der  Pol  term.  The 
reverse  statement  holds  for  the  other  solution,  the  amplitude  is  large  with  the  stabilization 
dictated  by  the  van  der  Pol  force  with  little  effect  of  friction.  As  the  coefficient  of  friction 
increases,  these  two  solutions  become  closer  together  and  eventually  merge.  The 
amplitude  of  response  then  appears  to  increase  monotonically  with  further  increases  of 
the  coefficient  of  friction.  The  effects  of  a  second,  viscous  damper  in  parallel  to  the 
friction  element  was  finally  considered.  It  was  found  that  the  two  damping  mechanisms 
do  not  necessarily  reinforce  each  other  and  that  the  friction  hinders  the  viscous 
dissipation  when  stick  occurs.  In  fact,  nonzero  amplitude  limit  cycle  oscillations  were 
noted  with  friction  in  situations  where  the  equilibrium  (zero  response)  would  in  fact  be 
stable  without  friction. 

Two  tightly  coupled  aeroelastic  systems  were  considered  in  task  (2): 

(A)  a  flat-plate  airfoil  model  supporting  a  torsional  friction  device  composed  of  a 
disk  flexibly  connected  to  the  plate  by  a  torsional  spring  and  squeezed  between  two  rough 
surfaces,  The  behavior  of  this  system  is  studied  when  placed  in  a  uniform,  inviscid  and 
incompressible  flow. 


(B)  a  NACA0012  airfoil  placed  in  a  uniform  inviscid  and  incompressible  air  flow 
and  supporting  either  the  same  frictional  device  as  the  flat  plate  or  a  block  sliding  in  a 
rough  internal  track. 

The  system  A  does  not  exhibit  limit  cycle  oscillations  in  the  absence  of  friction,  it 
is  either  stable  or  displays  an  explosive  flutter.  The  system  B  has  been  shown  in  some 
earlier  studies  to  exhibit,  in  the  absence  of  friction,  aerodynamic-driven  limit  cycle 
oscillations  near  the  flutter  speed.  The  analysis  of  these  two  systems  thus  provide  distinct 
perspectives  on  the  role  of  friction,  on  systems  exhibiting  explosive  flutter  (system  A) 
and  those  in  which  an  aero-driven  limit  cycle  oscillation  occurs  (system  B). 

The  system  A  was  analyzed  first  and  a  stability  analysis  of  it  was  initially 
conducted  with  the  disk  in  either  continuous  sticking  or  frictionless  slipping  modes  to 
assess  the  expected  stability  domains.  It  was  shown  and  justified  that  the  system  in  slip 
mode  exhibits  instabilities  at  earlier  flow  velocities  than  its  stuck  counterpart.  This 
property  allows  for  the  existence  of  both  super-  and  subcritical  limit  cycle  oscillations. 
While  the  subcritical  limit  cycles  were  observed  to  be  unstable,  a  zone  of  stable 
supercritical  limit  cycle  oscillations  was  found  that  extends  about  3%  past  the  flutter 
speed  of  the  system  without  the  friction  device.  This  gain  shows  a  good  stabilization 
property  since  the  moment  of  inertia  of  the  selected  friction  device  system  is  only  5%  of 
the  moment  of  inertia  of  the  plate.  The  observed  limit  cycles  exhibit  either  continuous 
slip  or  stick  slip  behaviors  and  are  either  single  frequency  (periodic)  or  multiple 
frequency  (aperiodic  or  chaotic)  with  the  latter  ones  appearing  primarily  at  the  highest 
flow  speeds  and  for  the  highest  frequencies  of  the  torsional  friction  device.  The  above 
results  were  obtained  by  time  marching  the  plate  equations  of  motions  with  a  rational 


approximation  of  the  Theodorsen  function  but  a  harmonic  balance  approach  was  also 
developed  that  led  to  very  good  approximations  of  the  single  frequency  continuous  slip 
limit  cycle  oscillations. 

The  system  B  was  considered  next  with  a  friction  component  modeled  by  a  block 
moving  in  an  internal  track.  While  limit  cycle  oscillations  were  observed,  it  was  also 
shown  that  the  block  could  become  stuck  at  a  position  far  from  its  original  one  and  thus 
would  create  a  change  of  inertia  sufficiently  large  to  stabilize  the  airfoil.  This  effect  does 
not  involve  any  dissipation  due  to  friction  and  is  thus  not  relevant  to  the  present  effort. 
Accordingly,  this  frictional  model  was  not  considered  further  and  was  replaced  by  its 
torsional  counterpart  (as  in  the  flat  plate  analysis)  which  does  not  suffer  from  the  same 
defect.  The  results  of  time  marching  computations  demonstrate  that  friction  can 
substantially  decrease  the  level  of  the  limit  cycle  oscillations,  especially  with  a  low 
coefficient  of  friction,  but  that  increases  in  the  response  are  also  possible,  depending  on 
the  selection  of  the  natural  frequency  of  the  torsional  friction  device.  As  in  the  flat  plate, 
continuous  slip  and  stick  slip  solutions  were  observed  most  of  which  were  single 
frequency  (periodic). 

The  last  task,  i.e.  (3),  of  this  project  focused  on  the  design,  fabrication,  and  testing 
of  an  actual  friction  device  on  an  airfoil.  The  basic  design  of  the  system  is  a  rotating  disk 
connected  to  the  airfoil  by  a  torsional  spring  of  natural  frequency  closely  matching  that  of 
the  airfoil  to  achieve  the  tuned  damper  arrangement  discussed  above.  The  spring  was 
instrumented  with  strain  gauges  to  measure  the  angle  of  torsion  of  the  system.  Finally,  a 
complex  system  was  designed  to  produce  the  required  normal  force  to  induce  friction  on 
the  moving  airfoil  without  transmitting  any  shear  or  moment.  Preliminary  shaker  testing 


in  the  ASU  vibrations  laboratory  successfully  validated  the  design:  stuck,  stick-slip,  and 
continuous  slip  motions  could  be  observed  by  varying  the  preload  on  the  disk.  Full  blown 
testing  in  the  DNW-TWG  transonic  wind  tunnel  at  DLR  led  to  the  observation  of  the  first 
recorded  limit  cycle  oscillations  with  friction  effects  and  relative  motions  of  the  disk  and 
airfoil.  Only  single  frequency  motions  were  observed  but  both  continuous  slip  and  stick 
slip  behaviors  were  found.  The  measurements  demonstrate  a  slight  effect  of  friction  and 
the  subcritical  nature  of  the  limit  cycle  oscillations.  This  behavior  was  explained  by 
analyzing  the  natural  frequencies  of  the  stuck  and  slipping  configurations  as  in  task  (2). 
Several  findings  are  however  still  unexplained,  e.g.  the  asymmetry  of  the  sticking  phase, 
the  continuous  slip  to  stick  slip  transition  associated  with  an  increase  in  response,  and  the 
near  Mach  independent  amplitude  of  torsional  response. 

The  overall  results  of  this  study  demonstrate  that  friction  can  indeed  provide  a 
stabilization  of  an  impending  flutter  and  can  significantly  decrease  the  amplitude  of 
existing  limit  cycle  oscillations  of  aeroelastic  systems  with  an  appropriate  selection  of  the 
friction  device  parameters  most  notably  natural  frequency  and  coefficient  of  friction. 
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ABSTRACT 


The  appearance  of  limit  cycle  oscillations  in  aeroelastic  systems  has  usually  been 
associated  with  nonlinearity  in  the  aerodynamics  and/or  in  the  structural  restoring  forces. 
It  has  however  been  recently  suggested  that  nonlinearity  in  the  damping  mechanism, 
more  notably  friction  between  a  small  moving  part  (or  many  such  parts)  and  the  wing, 
may  indeed  be  a  source  of  post-flutter  limit  cycle  oscillations.  The  present  work  provides 
a  first  validation  of  this  expectation. 

In  the  first  part  of  this  thesis,  a  two-degree-of-freedom  structural  dynamic  model 
is  considered  in  which  the  biggest  mass,  representing  the  wing,  and  the  smaller  one, 
modeling  the  sliding  component,  are  connected  by  at  least  a  friction  element.  The 
unstable  aerodynamic  effects  are  approximated  by  a  dashpot  of  negative  value  acting  on 
the  wing  alone.  The  response  of  this  system  is  analyzed  in  a  variety  of  situations,  i.e.  with 
an  additional  spring  between  masses,  an  additional  positive  damper  between  them,  and/or 
with  a  van  der  Pol  restoring  force  acting  on  the  wing  and  modeling  stabilizing  nonlinear 
aerodynamic  effects.  It  is  indeed  found  that  the  friction  can  stabilize  an  otherwise 
unstable  aeroelastic  system,  especially  when  the  natural  frequency  of  the  wing  and 
sliding  masses  considered  alone  are  close  together.  Periodic,  aperiodic,  and  chaotic 
motions  of  the  wing  are  observed  and  discussed.  The  presence  of  continuous  slip,  stick- 
slip,  and  stuck  motions  is  also  demonstrated  and  analyzed.  Finally,  predictions  from  the 
harmonic  balance  method  are  shown  to  match  well  the  amplitude  of  continuous  slip, 
periodic  motions. 

This  theoretical  work  was  complemented  by  a  series  of  tests  conducted  in  the 
transonic  wind  tunnel  of  the  Deutsches  Zentrum  fur  Luft  und  Raumfahrt  (DLR)  in 
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Gottingen,  Germany.  The  design  of  the  friction  test  article,  its  preliminary  validation  in 
the  ASU  vibrations  laboratory,  and  the  final  wind  tunnel  testing  are  presented  and 
discussed.  The  wind  tunnel  results  demonstrated  a  slight  stabilization  effect  of  friction. 
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CHAPTER  1 


INTRODUCTION 

Several  modern  aircraft,  most  notably  some  fighters  with  external  stores,  have 
been  found  to  exhibit  Limit  Cycle  Oscillations  (LCO),  i.e.  self-sustained  vibrations. 
Clearly,  these  vibrations  are  undesirable  from  all  perspectives:  they  affect  the  comfort  of 
the  pilot  and/or  passengers,  they  may  negatively  impact  the  accuracy  in  shooting,  and 
generally  produce  fatigue  in  the  aircraft  structure.  Accordingly,  a  series  of  investigations 
have  focused  on  explaining  the  occurrence  of  these  limit  cycle  oscillations  in  aircraft  and 
it  is  generally  recognized  that  these  vibrations  occur  in  the  neighborhood  of  flutter.  They 
further  have  been  classified  into  supercritical  (or  hard)  and  subcritical  (or  soft)  limit 
cycles  oscillations  with  the  critical  point  being  flutter.  Supercritical  LCO’s  thus  appear  at 
flow/flight  speeds  larger  than  the  flutter  speed  while  their  subcritical  counterparts  occur 
below  flutter  speed.  Both  supercritical  and  subcritical  LCO’s  are  the  manifestation  of  a 
nonlinear  aeroelastic  system,  as  a  linear  damped  aero-structure  model  can  only  exhibit  a 
decay  toward  rest  or  an  explosive,  divergent  behavior.  Thus,  paramount  to  LCO 
discussions  is  the  origin,  aerodynamic  and/or  structural,  and  nature  of  the  nonlinearity 
present. 

Nonlinear  aerodynamic  effects  have  often  been  invoked,  e.g.  the  transonic  shock 
separation  mechanism1,2  in  the  context  of  the  F-16,  but  structural  nonlinearities  have  also 
been  discussed.  Most  often  considered  have  been  the  large  deformation  stiffening  (as  in 
panel  flutter3),  and  freeplay4,5.  The  work  by  Chen  et  al.6  and  Mignolet  et  al.7  suggested 
however  that  a  different  mechanism,  specifically  friction,  could  qualitatively  explain 
many  of  the  LCO  findings  associated  with  the  F-16.  Interestingly,  such  an  explanation  is 
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also  supported  by  studies  performed  on  turbomachine  blades8,9  that  demonstrated  that 
friction  can,  in  certain  conditions,  stabilize  flutter  and  lead  to  supercritical  LCO.  A  recent 
investigation  has  also  shown  the  value  of  friction  to  damp  out  the  effects  of  gust10. 

In  this  light,  the  focus  of  the  present  investigation  and  of  a  related  one11  is  to 
conduct  an  in-depth  study  of  the  effects  of  friction  on  a  class  of  aeroelastic  systems 
exhibiting  either  explosive  flutter  or  limit  cycle  oscillations.  A  distinguishing  feature  of 
the  wing/wing-store  system  considered  here  is  that  the  friction  mechanism  must  be 
internal  to  the  system,  i.e.  no  rubbing  takes  place  between  components  separately 
supported  (e.g.  between  two  blades  of  a  rotor)  but  rather  the  components  in  relative  slip 
must  be  mounted  on  the  wing  (or  be  the  wing  itself).  In  keeping  with  the 
phenomenological  aspect  of  the  present  investigation,  the  analysis  was  conducted  under 
the  following  structural  simplifications. 

(i)  The  friction  effects  will  be  assumed  to  originate  from  a  single,  rigid 
component  moving  relatively  with  respect  to  the  wing  and  to  be  characterized  by  a 
preload  N  and  coefficients  of  friction  /us  and  fiD . 

(ii)  The  changes  of  external  structural  geometry  that  occur  due  to  the  relative 
sliding  will  be  neglected.  This  assumption  appears  justified  as  all  aircraft  components  are 
subjected  to  a  large  preload  so  that  a  large  amount  of  energy  can  be  dissipated  with  only 
small  relative  motions  at  the  joints.  Accordingly,  there  will  be  no  direct  effect  of  the 
moving  component  on  the  aerodynamics. 

Combining  the  above  two  assumptions  led  to  the  introduction11  of  an  airfoil 
model  with  internal  friction,  more  specifically  with  a  block  sliding  in  a  rough  track,  see 
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Fig.  1.1.  While  the  behavior  of  this  system  was  shown  to  exhibit  several  of  the 
characteristics  expected  from  friction,  it  was  also  found  to  be  very  sensitive  to  the  change 
of  inertia  that  occurs  as  the  block  is  moving  along  the  track.  In  fact,  it  was  noted  that  a 
static  displacement  of  the  block  in  the  track  provided  a  significant  stabilization  property 
which  rendered  extremely  difficult  the  assessment  of  the  effects  of  friction.  This  situation 
was  then  remedied12  by  assuming  that  the  rigid  component  of  Fig.  1.1  is  a  disk  that  spins 
around  its  center  of  mass  inside  the  airfoil,  see  Fig.  1.2.  Note  that  this  new  model  does 
satisfy,  as  well  as  its  counterpart  from  Fig.  1.1,  the  simplifications  (i)  and  (ii)  above. 


Fig.  1.1  Aeroelastic  system  with  sliding  block. 


Fig.  1.2  Aeroelastic  system  with  rotating  disk. 

In  this  overall  context,  the  specific  objectives  of  the  present  effort  are  twofold. 
First,  a  thorough  analysis  will  be  conducted  on  a  purely  structural  dynamic  model 
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representative  of  the  aeroelastic  systems  of  Fig.  1.2  to  assess  the  potential  of  friction  to 
stabilize  an  explosive  flutter  and  reduce  the  amplitude  of  existing  limit  cycle  oscillations. 
This  study  will  also  concentrate  on  the  properties  of  the  response  of  the  system  and  on  the 
influence  of  many  of  the  parameters  involved.  The  second  objective  is  the  design, 
fabrication,  and  testing  of  a  friction  device  similar  to  the  torsional  spring  and  disk  of  Fig. 
1.2.  The  preliminary  testing  of  the  device  was  successfully  accomplished  in  shaker  tests 
in  the  ASU  Vibrations  Laboratory.  Then,  its  final  testing  took  place  in  the  DLR  transonic 
wind  tunnel  mounted  on  a  NLR  7301  airfoil  undergoing  limit  cycle  oscillations  around 


Mach  number  0.6. 


CHAPTER  2 


TWO-DEGREE-OF-FREEDOM  DYNAMIC  MODEL 

The  focus  of  this  part  of  the  investigation  is  on  gaining  a  basic  physical 
understanding  of  the  effects  of  friction  on  the  stability/response  of  aeroelastic  systems. 
Given  that  very  little  is  known  about  this  topic  and  that  friction  inevitably  brings  out 
complex  nonlinear  characteristics,  it  was  decided  here  to  select  the  simplest  dynamical 
model  that  would  exhibit  the  features  of  the  problem.  Friction  in  aircraft  must  take  place 
between  two  components  of  the  aircraft,  as  opposed  to,  for  example,  the  wing  and  the 
ground.  A  two-degree-of-ffeedom  thus  seemed  the  lowest  order  system  to  be  considered. 
The  components  of  a  wing  are  well  fastened  together  and  thus  one  expects  that  a  relative 
motion  can  only  be  of  a  small  component.  Thus,  the  “sliding”  mass  ( M2  below)  must 
have  a  mass  much  less  than  the  wing  mass  (A/)).  It  remains  to  discuss  the  “aeroelastic” 
modeling.  In  keeping  with  the  novelty  of  the  topic,  it  was  decided  to  construct  a 
mathematically  simple  aerodynamic  model  that  would  exhibit  the  noted  explosive  flutter 
and  limit  cycle  oscillation  behavior.  These  possible  behaviors  are  included  in  the  model 
through  a  dashpot  of  negative  value  in  parallel  with  a  van  der  Pol  restoring  force.  The 
negative  dashpot  allows  unstable  motions  to  grow,  as  in  a  explosive  flutter  case,  but  the 
van  der  Pol  force,  if  present,  reinstates  stability  at  sufficiently  high  levels  of  response. 
Finally,  deformations  of  the  components  were  allowed  by  the  inclusion  of  springs.  The 
general  model  to  be  considered  here  is  thus  represented  in  Fig.  2.1. 

In  a  series  of  papers,  Tondl13'15  investigated  a  more  complex  version  of  Fig.  2.1  in 
which  a  structural  damper  also  exists  between  the  two  masses.  His  approach,  based  on 
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harmonic  balance  (see  section  2.4),  presented  only  limited  information  and  was  restricted 
by  the  necessary  assumptions  of  periodic  and  continuous  slip  motions.  Unfortunately,  as 
will  be  found  in  the  ensuing  sections,  these  assumptions  are  applicable  only  in  an 
extremely  limited  set  of  cases  so  that  the  results  of  References  11-13  neither  carry  the 
generality  nor  the  accuracy  required  here. 

It  should  finally  be  noted  that  the  structural  components  of  the  models  of  Figs  1.2 
and  2.1  (i.e.  without  aerodynamics  and  without  the  dashpot  c  and  van  der  Pol  term  5)  are 
in  fact  very  much  alike,  they  are  identical  if  the  plunging  motion  (vertical  displacement) 
of  the  airfoil  is  neglected.  In  this  case,  the  displacements  xj  and  X2  denote  the  pitching 
and  disk  angles,  the  masses  and  M2  are  in  fact  the  corresponding  inertias,  and  kx  is 
the  torsional  spring  stiffness. 


Fig.  2.1  2DOF  system  with  full  coupling 


2.1  Computational  Details 

The  system  in  Fig.  2.1  would  be  considered  a  linear  system  if  it  did  not  contain 
the  Coulomb  friction  element.  Friction  adds  significant  complexity  to  the  simulation  of 
the  system  in  that  it  can  completely  or  partially  inhibit  the  relative  motion  of  the  two 
masses  and  that  it  switches  direction  abruptly  as  the  relative  velocity  changes  sign.  In 
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fact,  during  unidirectional  (no  velocity  reversal)  slip  phases  as  well  as  stick  phases,  the 
equation  of  motion  of  the  system  are  linear,  so  that  the  nonlinearity  is  wholly 
concentrated  at  the  transitions.  These  observations  emphasize  the  need  to  accurately 
capture  the  transition  of  states  (stick-slip,  slip-stick,  slip-slip).  To  this  end,  the  time  step 
was  reduced  in  the  neighborhood  of  these  transitions  by  a  factor  as  large  as  128  through 
successive  halving.  Further,  the  equations  of  motion  were  integrated  numerically  using 
the  IMSL  routine  DIVPRK  (Runge-Kutta  of  orders  4-5)  with  a  low  tolerance  of  10'8 
(both  relative  and  absolute  errors)  and  a  time  step  typically  smaller  or  equal  to  0.01  s 
while  the  fundamental  linear  frequency  is  close  to  1  radian  per  second. 

2.2  Assessment  of  Friction  on  an  Aeroelastic  System  Exhibiting  Explosive  Flutter 

The  analysis  of  the  2-degree-of-freedom  system  was  conducted  in  two  separate 
steps.  It  was  first  analyzed  under  the  assumption  5  =  0  which  implies  that  no  nonlinear 
stabilizing  aerodynamic  effect  exists.  Accordingly,  the  linear  system  (without  friction) 
exhibits  only  an  unstable  behavior  typical  of  explosive  flutter.  The  considerations  of 
stabilizing  nonlinear  aerodynamics,  i.e.  5*0,  will  be  conducted  in  Section  2.3.  Two 
separate  cases  were  further  considered  with  5  =  0,  i.e.  k{=0  or  k,^0.  Under  the  first 
condition,  the  2-degree-of-freedom  system  obtained  in  slip  phases  is  classically  damped 
(at  the  contrary  of  the  situation  kx* 0)  and  an  exact  formulation  (EF)  of  the  response  can 
be  derived  without  too  much  difficulty.  This  approach,  as  well  as  a  direct  numerical 
simulation,  (NS),  of  the  equations  of  motion  were  thus  performed  for  kx=0 .  In  the  second 
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step,  nonzero  values  of  kx  were  considered  but  only  the  numerical  simulation  was  carried 
out. 

2.2.1  Primary-Secondary  Coupling  Through  the  Friction  Element  Only 

The  system  under  consideration  in  this  section  is  the  2-degree-of-freedom  (2DOF) 
system  shown  in  Fig.  2.2. 


Fig.  2.2  2DOF  system  with  friction  coupling. 


2.2. 1.1  Exact  Formulation 

Following  Den  Hartog16  (see  also  Pesheck  and  Pierre17),  a  periodic  solution  of  the 
equations  of  motion  can  be  obtained  by  marching  the  solution  from  unknown  initial 
conditions  at  t  =  0  to  an  unknown  time  T  at  which  periodicity  conditions  are  enforced. 
Such  an  approach  can  be  used  for  continuous  slip  behavior  as  well  as  solutions  with 
stick-slip  but  the  number  of  stick-slip  phases  must  be  chosen  a  priori.  This  process 
transforms  the  solution  of  the  nonlinear  differential  equations  into  the  solution  of  a  set  of 
nonlinear  algebraic  equations  for  the  unknown  initial  conditions,  period  (or  half  period 
T ),  and  stuck  times  t,  if  applicable.  This  process  is  exemplified  below  on  the  system  of 


Fig.  2.2. 
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Since  the  equations  of  motion  of  the  system  of  Fig.  2.2  are  time  invariant,  the  time 
t  =  0  can  be  arbitrarily  selected.  Then,  assume  for  convenience  that  this  time  corresponds 
to  a  stick-slip  or  slip-slip  transition.  Thus,  the  initial  conditions  were  selected  as 

*i  (°)  =  *2  (°)  =  v0 ;  *i  (0)  =  x10  and  x2  (o)  =  x20 .  (2-1) 

During  the  ensuing  slip  phase,  it  is  assumed  that  x,>x2  so  that  the  response  is  described 
by  the  equations  of  motion 

Mxxx  +  cx,  +  fct,  =  Fix  and  M2x2  =  -Fn  (2-2) 

where 


F31  =  JJdN  sgn(x2  -x,) 

The  response  of  the  two  blocks  can  thus  be  expressed  as 

xl(t)  =  e<i“'1'[A2  cos a>’dlt  +  B2  sin co'd2t] - 

k 


and 


where 


(/)  =  ■i— •  f 2  +  v0t  +  x. 
w  2  M,  0 


^  ~  2^'  ;  dd2  ^ 


,  __  ,  MdN.  n  _  v0  +O2A2 

a2~xiq  +  ~ r~>  v-l  - ; - 

k  co 


d  2 


(2-3) 


(2-4) 


(2-5) 


(2-6) 


(2-7) 


If  a  continuous  slip  solution  takes  place,  the  unknown  initial  conditions  in  Eq.  (2-1)  and 
the  half-cycle  time  T  will  be  evaluated  from  the  periodicity  conditions 

xi{T)=-x\{°)  =  -xw'’  x,(r)  =  -Xi(o)  =  -Vo;  j2(r)  =  -x2(0)  =  -v0. 


(2-8) 
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Enforcing  these  constraints  leads  to  a  nonlinear  equation  for  the  half-cycle  time  T  and 
ensuing  linear  equations  for  the  unknown  initial  conditions.  Specifically,  imposing  the 
periodicity  condition  at  T  for  x2  yields 

~~~  T  =  -2v0 .  (2-9) 

M2 

The  periodicity  conditions  for  x,  and  jc, ,  respectively,  yield 


A,  sin co'j-, T  =v„  e^^  +  cos co\,T C'— sin a>'nT 


(2-10) 


A2  e~CiaiT  cos eo'dlT  +  -’2-(°^-svixco'd2T  +1  +  -^-e  (l<0lT  sinter  =  -  .(2-11) 

L  V  <°d 2  J  \  <2  k 

Combining  Eq.  2-9-2-11  then  results  in  the  single  nonlinear  algebraic  equation  for  the 


half-period  T 


2r  sin  0)'d2T  +  co'd2T[coscod2T  +  cosh£2ci)2T]  =  0. 


(2-12) 


A  stick-slip  solution  was  also  obtained.  The  transition  from  slip  to  stick  was 
assumed  to  take  place  at  time  t  =  t0  such  that  the  relative  velocity  of  the  two  blocks 


vanishes,  i.e. 


*i(0- *2(0=0 


(2-13) 


[v0  cos o)dlt0  -{C2(o'2B2  +  A2a)d2 )sin a)d2tQ ]  =  10  +  v0 .  (2-14) 


After  the  transition,  the  slider  sticks  and  the  system  is  governed  by  the  equations  of 


motion 
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(M,  +  M2)x1  +  ci,  +  fee,  =  0  and  ft,  =  x2 
The  corresponding  response  can  thus  be  expressed  as 

x,  (t)  =  e~Cia'‘  [ax  cos  co’dlt  +  Bl  sin  co'dlt\ 

where 


(2-15) 

(2-16) 

(2-17) 


2.2. 1.2  Numerical  Results 

A  parametric  study  of  the  2-degree-of- freedom  system  with  ft,  =0  has  revealed 
that  the  internal  friction  mechanism  stabilizes  the  unstable  aerodynamics  and  leads  to 
LCO  for  a  broad  range  of  negative  damping  ratios  (see  Fig.  2.3).  Further,  the  response  of 
the  system  was  found  to  exhibit  many  of  the  features  already  observed  in  connection  with 
the  1.5  degree-of-freedom  system6'7,17,  e.g.  at  most  one  stable  solution  was  observed  in  all 
cases  considered,  an  upper-lower  branch  structure  exists  in  the  plot  of  amplitude  versus 
damping,  Fig.  2.4  (the  lower  branch  is  stable  and  denoted  by  S  while  the  upper  branch  is 
unstable  and  denoted  by  U),  and  the  solutions  were  both  periodic  and  symmetric. 
However,  at  the  contrary  of  that  simpler  system,  the  two-degree-of-freedom  model  with 
ft,  =  0  displays  two  types  of  LCO  solutions:  continuous  slip  (CS)  and  stick-slip  (SS,  a 
single  stick  phase  per  half-cycle). 
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While  the  unstable  solutions  were  all  characterized  by  a  continuous  slip  behavior, 
the  stable  motions  were  primarily  corresponding  to  a  stick-slip  behavior  although  a  small 
zone  of  stable  continuous  slip  responses  were  observed  near  the  largest  stabilizable 
damping  ratio.  The  results  obtained  with  the  exact  formulation  matched  very  closely  their 
numerical  simulation  counterparts  even  in  a  stick-slip  motion18. 

The  results  shown  in  Fig.  2.3  and  2.4  were  obtained  for  a  typical  friction  setup 
(/i5=uD=0.25 ,  iV  =  l,  k-  1)  but  for  a  rather  large  mass  ratio,  i.e.  Mx=\ ,  M2- 0.25. 
While  the  phenomenological  features  described  above  were  found  to  also  hold  at  much 
lower  values  of  the  secondary/sliding  mass  M2 ,  it  was  found  that  the  largest  stabilizable 
damping  ratio  is  approximately  proportional  to  the  mass  ratio,  r-M2/Ml.  A  relationship 
between  the  damping  and  mass  ratios  is  now  desired  to  fully  understand  Fig.  2.4.  To  this 
end,  the  condition  of  Eq.  (2-12)  was  reconsidered  and  it  was  argued  that  the  largest 
damping  ratio  is  such  that  (see  also  Fig.  2.4) 

—  =  oo  or  =  (2-18) 

d#  dr 

Then,  differentiating  Eq.  (2-12)  with  respect  to  T  and  using  Eq.  (2.18)  leads  to 

2r  cos  co'dlT  +  cos  co'd2T  +  cosh  {£'2co'2T) 

+  co'd2T  [-  sin  a)d2T  +  sinh {C2o)'d2T )]  =  0. 

Equations  (2-12)  and  (2-19)  can  be  simplified  as 

2r  sin  u  +  w[cosw  +  cosh  //h]  =  0 

2r  cos  u  +  cos  u  +  cosh  rju-u  sin  u  +  tju  sinh  jju  =  0 


(2-19) 

(2-20) 

(2-21) 


where 
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=  jju  .  (2-22) 

It  is  then  desired  to  solve  Eq.  (2-20)  and  (2-21)  for  the  damping  ratio  and  half-period,  i.e. 
for  u  and  77  given  r  small.  When  the  mass  ratio  goes  to  zero,  it  is  expected  that  the 
maximum  stabilizable  damping  is  zero  and  thus  that  there  exists  a  solution  of  Eq.  (2-20) 
and  (2-21)  with  r  -  77  =  0.  It  is  indeed  the  case  and  requires  that  u  =  n.  For  small  but 
nonzero  mass  ratios,  a  perturbation  analysis  will  be  carried  out  with 

u  =  n  +  a .  (2-23) 


U~®d2?  >  £2^2^  ~ 


a 


In  fact,  it  is  found  that 


-2  ra  +  ^-[a2  +tj2x2)=  0, 

(2-24) 

2r 

a-  — . 

(2-25) 

7C 

In  terms  of  the  half-period  and  damping  ratio,  these  results  imply 

T  ±2r 

T  =  7T  + - 

K 

(2-26) 

and 


2  JkM\  J 


2  r 


7C 


for  r«l . 


(2-27) 


Half  Period,  T  Amplitude 
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According  to  this  relation,  a  2%  damping  ratio  stabilization  would  require  a  mass 
ratio  r  of  at  least  0.1,  a  larger  value  than  would  be  expected  to  be  typical  of  the  F-16.  A 
larger  stabilization  capacity  of  the  two  mass  system  of  Fig.  2.2,  and  thus  potentially  a 
lower  mass  ratio,  would  require  an  increased  relative  response  of  the  two  masses.  Such  an 
increase  can  be  achieved  by  increasing  the  transfer  of  energy  between  the  two  masses.  In 
view  of  the  primary-secondary  nature  of  the  system,  it  is  expected  that  the  difference  in 

the  natural  frequencies  of  the  primary  (a>]  =  )  and  the  secondary  (co2= 0  for  the 

system  with  Ar,=0)  will  be  the  dominant  parameter  controlling  the  energy  transfer. 
Further,  a  maximum  transmission  could  be  conjectured  to  occur  when  cox&co2 . 

2.2.2.  Primary-Secondary  Coupling  Through  the  Friction  Element  and  a  Spring 


Fig.  2.5  2-degree-of-freedom  system  with  friction  and  stiffness  coupling. 

In  this  section,  the  stiffness  k{  was  varied  from  zero  and  the  system  response  was 
obtained  by  numerical  integration.  A  series  of  values  of  the  stiffness  kx  and  of  the 
damping  coefficient  c  were  considered  with  the  mass  ratio  r  =  0.05  and  the  single  set  of 
initial  conditions  Xx  (o)  =X2  (o)  =  0. 1  and  Xx  (o)  =X2  (o)  =  0 . 
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The  wealth  of  results  obtained  (see  Fig.  2.6)  has  shown,  first  and  foremost,  that 
the  above  conjectures  were  perfectly  valid.  Indeed,  each  of  the  dots  on  Fig.  2.6(a) 
indicates  a  stable  solution  and  thus: 

1)  an  increased  stabilization  capability  is  possible  by  varying  the  natural  frequency  of 
the  sliding  mass  alone,  and 

2)  the  maximum  stabilization  possible  appears  to  take  place  when  co^-co^  or 
equivalently  k{=r=  0.05  here  since  k  is  set  to  unity. 

In  analyzing  these  results,  it  could  be  argued  that  the  increase  in  stabilization 
could  simply  be  a  result  of  an  increase  in  stiffness  and  a  corresponding  reduction  of  the 
(negative)  damping  ratios  of  the  system  since  it  is  the  dashpot  value  c  which  is  fixed. 
These  comments  suggest  that  a  more  definite  proof  of  the  increased  stabilization  would 
require  the  analysis  of  the  most  negative  damping  ratio  for  each  of  the  stable  systems  of 
Fig.  2.6(a).  This  effort  yielded  Fig.  2.6(b),  which  also  displays  the  trends  stated  above, 
thereby  demonstrating  their  validity.  Note  further  that  the  largest  stabilized  damping  ratio 
observed  was  3.5%  which  is  almost  twice  what  is  required  for  the  F-166  while  the  mass 
ratio  of  0.05  is  quite  reasonable.  It  is  then  concluded  that  the  proposed  friction  based 
LCO  mechanism  is  indeed  possible. 

Having  established  the  existence  of  an  increased  domain  of  stable  solutions,  it  is 
important  to  assess  the  phenomenological  characteristics  of  these  solutions  and  to 
compare  them  with  those  obtained  with  kx= 0 .  In  fact,  four  different  types  of  solutions 
were  observed  when  k^p 0  that  can  be  characterized  as  single  frequency  stick-slip  (SFSS), 
single  frequency  continuous  slip  (SFCS),  multiple  frequency  stick-slip  (MFSS),  and 
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multiple  frequency  continuous  slip  (MFCS).  The  term  single  frequency  is  used  here 
somewhat  freely  to  describe  a  solution  with  a  simple  period  and  with  minimal  harmonics, 
see  Fig.  2.7  for  such  a  solution  (continuous  slip  example).  Note  in  particular  that  the 
magnitude  of  the  third  harmonic  of  the  response  of  the  mass  M,  is  only  1%  of  its 
fundamental  counterpart.  For  comparison,  “multiple”  frequency  solutions,  see  Figs  2.8 
and  2.9,  exhibit  a  beating  phenomenon  (see  Figs  2.8(a)  and  2.9(a))  that  is  associated  with 
the  presence  of  two  frequencies  close  to  the  dominant  one,  see  Fig.  2.8(b).  The 
corresponding  phase  plane  plots  (Figs  2.8(b),  2.9(c))  clearly  show  the  complexity  of  these 
motions  as  compared  to  the  single  frequency  solution,  see  Fig.  2.7(c).  A  different 
perspective  into  the  phase  plane  features  can  be  obtained  by  sampling  the  time  histories  at 
specific  times.  In  forced  response  problems,  the  sampling  is  accomplished  every  cycle  of 
the  excitation  but  in  the  present  self-excited  problem,  it  was  found  convenient  to  monitor 
the  displacement  and  velocity  of  mass  Mx  when  the  relative  motion  achieved  maximum 
values.  The  corresponding  sampled  phase  plane  plots  are  shown  in  Fig.  2.10.  For  a 
periodic  motion,  only  a  single  dot  appears,  as  seen  for  the  SFCS  solution.  The  MFCS 
solution  on  the  contrary  leads  to  a  close  curve  while  the  MFSS  solution  is  characterized 
by  a  similar  curve  and  a  densely  populated  domain.  It  is  tentatively  concluded  from  those 
figures  that  the  corresponding  MFCS  and  MFSS  solutions  are  aperiodic  and  chaotic, 
respectively.  Note  that  the  difference  in  behavior  between  kx  —  0.039  and  0.040  was  not 
easily  detectable  from  the  time  histories  and/or  full  phase  plane  plots  but  is  quite  clear 
from  the  Poincare-type  plot  of  Fig.  2.10. 


Dashpot  constant 
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It  was  next  desired  to  analyze  the  “magnitude”  of  the  response  as  a  function  of  the 
stiffness  kx  and  the  damping  constant  c.  Since  an  amplitude  cannot  be  properly  defined 
in  the  multiple  frequency  situation,  it  was  decided  to  study  the  highest  level  of  response 
obtained  in  the  steady  state  part  of  the  computed  response  instead,  see  Fig.  2.11.  Note 
that  the  curves  are  not  plotted  continuously  at  the  change  of  solution  type,  i.e.  from  single 
to  multiple  frequency  and  vice  versa.  The  general  tendency  of  these  curves  is  that  the 
level  of  response  decreases  from  the  edges  of  the  stable  zone  toward  its  center 
kx  «  r  =  0.05 ,  as  expected  from  the  increase  in  energy  transfer.  Further,  the  largest 
response  is  substantially  increased  by  the  beating  phenomenon  of  Figs  2.8(a)  and  2.9(a), 
i.e.  as  the  solution  switches  from  a  single  to  two  dominant  frequencies. 

2.2.3  Changes  in  Solution  Types  -  Bifurcations 

The  arrangement  of  the  4  solutions  types  in  the  c  -  kx  or  -kx  domain,  see  Fig. 
2.6,  is  very  interesting.  Near  the  edges  of  the  stable  zone,  it  appears  that  the  single 
frequency  solution  only  is  present.  Further,  it  exhibits  a  stick-slip  (SFSS)  behavior  only 
for  the  smallest  values  of  the  negative  damping  and  a  continuous  slip  (SFSS)  through 
most  of  the  domain.  This  observation  is  physically  expected  as  a  continued  slip  is  likely 
to  provide  an  increased  dissipation  and  in  that  respect  is  consistent  with  the  kx  =  0 
findings  of  Fig.  2.4(a).  As  the  stiffness  kx  is  varied  (increased  or  decreased)  toward  the 
center  of  the  stable  domain  (kx»r),  a  change  of  solution  takes  place  and  the  single 
frequency  solution  switches  to  a  multiple  frequency  one,  first  maintaining  the  continuous 
slip  (for  sufficiently  large  negative  damping),  then  exhibiting  stick  phases. 
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The  appearance  of  the  multiple  frequency  solutions,  not  seen  in  the  simpler, 
=  0  case,  may  be  associated  with  the  two  natural  frequencies  of  the  linear  system. 
Indeed,  for  small  mass  ratios  and  for  values  of  kx  ~  r,  the  two  system  frequencies  are 
close  to  each  other  and  could  be  producing  the  tightly  spaced  peaks  of  Fig.  2.9(b). 
Further,  when  kx  is  exactly  equal  to  r,  the  corresponding  (negative  here)  damping  ratios 
are  also  equal  so  that  contributions  of  both  linear  modes  to  the  response  would  be 
certainly  be  expected  in  this  case.  As  the  stiffness  kx  deviates  from  the  mass  ratio  r,  the 
two  damping  ratios  rapidly  differ  from  each  other  so  that  the  response  should  emphasize 
more  heavily  one  of  the  two  modes,  i.e.  the  one  characterized  by  the  most  negative 
damping  ratio.  These  observations  provide  a  tentative  explanation  for  the  presence  of 
three  dominant  frequencies  near  the  center  of  the  stable  domain,  i.e.  «  r,  but  only  one 
frequency  near  the  edges  of  this  region  as  seen  in  Fig.  2.6. 


a)  Normalized  Fourier  Transform  of  Xx 

Fig.  2.7  A  typical  "single"  frequency  continuous  slip  solution,  M1  =  0.05, 

kx  =  0.06,  c  =  -0.07. 


Velocity 
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Fig.  2.11  Largest  response  of  mass  Mx  as  a  function  of  the  stiffness  kx  for  different 

values  of  the  damping  constant  c. 

A  more  formal  approach  can  be  undertaken  to  better  clarify  the  transition, 
bifurcation  in  fact,  of  solutions  from  single  to  multiple  frequencies.  This  approach, 
referred  to  as  Floquet  analysis19,  relies  on  the  availability  of  the  periodic,  steady  state 
solution  at  a  time  t,  x(/).  If  the  equations  of  motion  can  be  linearized  around  this  steady 
state  solution,  the  response  after  one  cycle,  x(f)  +  8x(/  +  T),  corresponding  to  an  imposed 
perturbations  at  the  beginning  of  the  cycle,  i.e.  5x(i),  can  be  expressed  as 

5x{t  +  T)  =  P5x{t).  (2-28) 

In  this  equation,  P  is  a  square  transition  matrix  which  can  be  identified  by  proceeding  in 
turn  with  small  perturbations  of  each  of  the  state  variables  (4  in  the  present  context)  and 
computing  the  corresponding  perturbed  solutions  after  one  cycle.  Repeating  the  argument 
of  Eq.  (2-28),  the  perturbation  after  n  cycles  is  related  to  the  original  one,  5x(/),  by  the 

matrix  Pn .  Then,  the  steady  state  solution  x(t)  will  be  stable,  i.e.  the  effects  of  the 
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perturbation  8x(t )  will  disappear  as  t  — >  oo  ,  if  Pn  ->  0  when  n  — >  co  or  equivalently  if 
the  eigenvalues  of  P  are  all  less  than  1  in  magnitude. 

This  result  is  valid  in  forced  vibrations  cases  in  which  time  is  defined  by  the 
excitation.  In  self  excited  problems  however,  it  is  possible  for  the  perturbed  solution  to 
converge  to  the  unperturbed  steady  state  response  but  time  shifted  by  8/ .  In  fact,  this 
situation  occurs  for  8x(/)  =  x(t)  8t  as  *(/)  +  8 x(t)  is  then  simply  x(t  +  8/) .  Thus,  after 
one  period,  the  perturbed  response  is  x(t  +  8t  +  T)  =  x(t  +  8 1)  from  the  periodicity  of  the 
exact  steady  state  solution.  From  this  relation,  it  is  thus  concluded  that  an  eigenvalue  of 
+1  should  always  be  present  in  P  in  self  excited  problems. 

The  estimation  of  the  matrix  P  was  performed  for  c  =  -0.07  and  a  series  of  values 
of  k j  for  which  a  single  frequency  solution  is  encountered,  i.e.  k\  e  (0.34,  0.3885]  and  the 
corresponding  eigenvalues  (complex  and  imaginary  parts)  are  plotted  in  Fig.  2.12.  The 
persistent  eigenvalue  at  +1  is  clearly  seen  and  so  is  the  convergence  of  two  eigenvalues 
toward  the  unit  circle  as  k\  ->  0.3885.  From  this  last  observation,  it  is  tentatively 
suggested  that  the  transition  from  a  single  to  multiple  frequencies  is  a  Hopf  bifurcation. 

Considering  this  data  for  monotonically  decreasing  values  of  k\  provides  also  a 
perspective  on  the  loss  of  stability  that  occurs  for  k\  ~  0.34 .  Specifically,  it  is  noted  that 
the  real  eigenvalue  of  value  <1  is  slowly  progressing  toward  the  unit  circle,  suggesting 
that  the  loss  of  stability  is  accomplished  through  the  occurrence  of  a  double  eigenvalue  at 
+1. 
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Fig.  2.12  Eigenvalue  plot  of  Floquet  analysis,  c  =  -0.07. 

2.3  Assessment  of  Friction  on  an  Aeroelastic  System  Exhibiting  Post  Flutter 
LCO 

2.3.1  Primary-Secondary  Coupling  Through  the  Friction  Element  and  a  Spring 


Fig.  2.13  2-degree-of-freedom  system  with  van  der  Pol  force. 
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The  simple  system  of  Fig.  2.5  does  not  exhibit  limit  cycle  oscillations  in  the 
absence  of  friction,  its  response  exponentially  increases  as  time  grows.  It  was  thus 
necessary  to  add  another  stabilizing  factor  to  induce  limit  cycle  oscillations  in  the 
absence  of  friction.  One  such  term  often  encountered  is  a  van  der  Pol  (vdP)  force  which 
is  proportional  to  the  product  of  the  square  of  the  displacement  and  the  velocity.  This 
additional  force  was  assumed  to  act  on  the  primary  mass,  see  Fig.  2.13,  yielding  the 
equations  of  motion 

Mx  jC[  +c(l-£ xx  +  (k+kx  )x,  -kx  x2  =  F2X  (2-28) 

and 

M2x2  -kx  x,  +kx  x2  =  -F31  (2-29) 


where 


Fn  =  judN  sgn(x2  -Xj)  during  slip 
|FjJ  <  /jsN  during  stick. 


(2-3  Oa) 
(2-30b) 


A  nondimensionalization  of  the  equations  of  motion  provides  significant  insight  into  the 
behavior  of  the  system.  To  this  end,  introduce 

cox=^Jk/ Mx  ;  co2=^jkx  / M2  ;  r-M2  / Mx ; 

T=coxt ;  yx=xx(k  /  fiD N) ,  i  =  1 ,2; 

C=ccox  /  2k  ;  'S=5k1  /  fi^N2 ;  jfe,  =*,/*.  (2-31) 

Then,  the  equations  of  motion,  Eq.  (2-28)-(2-30),  can  be  rewritten  as 

y:+2c{l-Sy?)y'x  +{l+ki)yx  -hy2  =F*  (2-32) 


and 
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ry2~kiyl+kiy2=-F2\ 

where  (.)'  denotes  the  derivative  with  respect  to  t, 


and 


Fi\  =  sgn(j4  -  _y{)  during  slip 


F  31  <  /js  I  nD  during  stick. 


(2-33) 


(2-34a) 


(2-34b) 


Note  in  particular  from  Eq.  (2-32)-(2-34)  that  the  amplitude  of  response  is 
proportional  to  juD  if  /js  -/jd  and  S  =  0 .  Further,  this  property  also  holds  for  nonzero 


values  of  <5  provided  that  8  is  held  constant. 

As  in  the  previous  section,  a  numerical  approach  was  adopted  to  study  the 
response  of  the  system.  As  before,  Mx  and  M2  were  selected  to  be  1  and  0.05,  and  k  was 
set  to  1.  Further,  the  static  and  dynamic  coefficients  of  friction  were  selected  equal, 
Ms  ~Md  -  M-  The  other  parameters  were  varied  to  provide  a  thorough  perspective  on  the 
response  of  the  system. 

The  system  of  Fig.  2.13  was  analyzed  with  different  values  of  the  van  der  Pol 
constant  8  and  the  amplitude  of  response  and  stability  plots  for  each  of  these  values  was 
generated.  Shown  in  Fig.  2.14-2.17  are  the  results  corresponding  to  <5=0.025  (Fig.  2.14 
and  2.15)  and  <5=0.04  (Fig.  2.16  and  2.17).  To  clarify  the  discussion  of  these  figures, 
note  first  that  the  van  der  Pol  system  is  stable  for  all  cases  considered  and  thus 
stabilization  must  occur  for  <5>  0 .  Further,  the  response  of  the  van  der  Pol  system  without 
friction  is  a  single  frequency  (continuous  slip)  solution. 
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Fig.  2.14  Largest  response  of  mass  Mx  as  a  function  of  the  stiffness  kx  for  different 
values  of  the  damping  constant  c,  8= . 025,  and  ju=  0.25. 


Fig.  2.15  Domain  and  type  of  stable  solutions  8  =  0.025  and  //=  0.25. 
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—  0.02  —  0.03  —  0.04  —  0.05  —  0.06  —  0.07  —  0.08 

—  0.02/2  —  0.03/2  —  0.04/2  0.05/2  —  0.06/2  —  0.07/2  —  0.08/2 


Fig.  2.18  Largest  response  of  mass  M,  as  a  function  of  the  coefficient  of  friction  m, 
for  different  stiffnesses  ft, ,  c  =  -0.04,  and  8  -  0.06  (The  notations  /2  refers  to  a 

second  solution  for  that  ft, ). 


—  0.02  : — 0.03  — 0.04  — 0.05  — 0.06  — 0.07  — 0.08 

—  0.02/2  — 0.03/2  — 0.04/2  —0.05/2  —0.06/2  —0.07/2  —0.08/2 


Fig.  2.19  Largest  response  of  mass  Mx  as  a  function  of  the  coefficient  of  friction  m, 
for  different  stiffnesses  ft, ,  c  =  -0.09,  and  5  =  0.06  (The  notations  /2  refers  to  a 

second  solution  for  that  ft, ). 
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Comparing  Fig.  2.6(a),  2.15,  and  2.17,  it  is  first  seen  that  the  unstable  regions  of 
Fig.  2.6(a)  are  replaced  by  stable  single  frequency  continuous  slip  motions  characteristic 
of  a  van  der  Pol  dominated  response.  Note  further  that  the  domain  of  multiple  frequency 
solutions  shrinks  as  5  is  increased  leaving  mostly  single  frequency  motions.  This  finding 
may  be  explained  by  noting  that  the  multiple  frequency  solutions  typically  lead  to  large 
amplitudes  because  of  the  beating,  see  Fig.  2.11.  For  such  large  amplitudes  however,  the 
stabilization  of  the  van  der  Pol  term  is  significant  forcing  a  decrease,  and  eventually 
disappearance,  of  the  beating  effect.  The  decreased  presence  of  the  multiple  frequency 
solutions  is  also  visible  from  the  amplitude  of  response  plots,  Fig.  2.11,  2.14,  and  2.16, 
where  they  imply  higher  amplitudes  than  the  neighboring  single  frequency  solutions. 
Notice  on  these  same  figures  that  the  sharp  transitions  to  instability  in  Fig.  2.11  are 
replaced  by  smooth  increases  toward  the  amplitude  corresponding  to  the  van  der  Pol 
alone. 

Another  change  in  the  character  of  the  solution  concerns  the  transition  from  stick- 
slip  to  continuous  slip.  It  can  be  seen  from  Fig.  2.6(a),  2.15,  and  2.17  that  stick-slip 
solutions  typically  occur  at  low  values  of  the  negative  damping  ratio  and  correspond  to 
smaller  response  level  for  which  the  dissipation  during  the  slip  phase  is  sufficient  to 
stabilize  the  motions.  Continuous  slip  solutions  on  the  contrary  appear  at  higher  response 
levels  where  the  amplitude  of  response  is  large  enough  to  force  the  continuous  motion  of 
the  sliding  block.  It  is  seen  from  Fig.  2.15  and  2.17  that  the  value  of  the  negative 
damping  ratio  at  which  the  stick-slip  to  continuous  slip  transition  occurs  increases  (in 
magnitude)  as  the  van  der  Pol  parameter  5  is  increased.  This  result  is  in  fact  expected  as 
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an  increase  in  5  promotes  stability  and  forces  lower  response  levels  thereby  favoring 
stick-slip  solutions.  An  increase  of  the  damping  constant  (larger  negative)  must  then  take 
place  to  increase  the  energy  introduced  in  the  system  and  to  induce  a  continuous  slip 
behavior. 

The  V-shape  of  the  amplitude  plot,  slightly  visible  for  <5=0  (Fig.  2.11),  quite 
clear  for  <5=0.025  (Fig.  2.14),  and  obvious  for  <5=0.04  (Fig.  2.16),  demonstrates  that  an 
amplitude  of  response  substantially  below  the  van  der  Pol  alone  level  can  be  achieved 
with  friction.  Thus,  friction  can  also  play  an  important  role  in  reducing  the  level  of 
response  of  existing  limit  cycle  oscillations.  Note  further  that  friction  is  most  effective 
when  &,«0.05,  especially  for  large  negative  damping  ratios,  as  was  already  noted  in  the 
previous  section. 

The  above  discussion  has  focused  on  the  behavior  of  the  response  as  a  function  of 
<5,  c,  and  k{  for  a  fixed  value  of  pi  (=  0.25)  and  it  is  now  desired  to  assess  the  role  of  the 
friction  coefficient.  To  this  end,  the  value  of  <5  was  set  to  0.06  and  the  response  of  the 
system  of  Fig.  2.13  was  determined  for  two  damping  coefficients,  i.e.  c  =  -0.04  and  - 
0.09,  as  function  of  k{  and  pi.  Surprisingly,  it  was  found  for  both  values  of  c  that  two 
stable  solutions  exist  for  low  coefficients  of  friction,  see  Fig.  2.18  and  2.19.  Each  of  these 
two  solutions  essentially  involves  a  single  stabilization  factor:  the  very  low  amplitude 
solution  is  associated  with  the  low  friction  coefficient  without  any  significant  effect  of 
the  van  der  Pol  term.  In  such  conditions,  a  stable  solution  must  exhibit  a  small  amplitude 
as  the  low  level  of  friction  will  not  provide  much  dissipation.  At  the  contrary,  the  large 
amplitude  solution  is  only  marginally  affected  by  the  small  frictional  term  and  it  is 
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primarily  the  van  der  Pol  force  that  dictates  the  response.  As  the  coefficient  of  friction 
increases,  its  effect  on  the  high  amplitude  solution  increases  implying  a  decrease  of  the 
response  level.  Further,  a  higher  coefficient  of  friction  provides  a  balance  of  the  unstable 
aerodynamic  forces  at  higher  amplitude  levels.  Thus,  the  amplitude  of  the  lower  (resp. 
higher)  solution  increases  (resp.  decreases)  until  they  both  merge  into  a  single  one 
solution  as  seen  in  Fig.  2.18  and  2.19.  The  value  of  the  coefficient  of  friction  at  which  the 
merging  of  the  solutions  occurs  is  a  particular  good  operating  point  as  it  yields  the 
smallest  value  of  the  maximum  stable  solution.  Note  further  that  the  existence  of  two 
solutions  occurs  at  high  negative  aerodynamic  damping  only  when  the  natural 
frequencies  of  the  two  blocks  alone  are  identical  or  almost  so. 

2.3.2  Primary-Secondary  Coupling  Through  the  Friction  Element,  a  Spring,  and  a 
Dashpot 

Dissipation  associated  with  the  relative  motion  of  the  two  masses  may  include 
more  than  just  friction,  most  notably  viscous  effects,  and  it  is  finally  desired  here  to 
assess  how  the  presence  of  an  additional  damping  mechanism  would  affect  the 
conclusions  drawn  in  connection  with  the  system  of  Fig.  2.13.  To  this  end,  a  dashpot 
(viscous  element)  was  added  between  the  two  masses,  see  Fig.  2.1,  and  a  parametric 
study  of  the  response  magnitude  and  character  (i.e.  single  vs.  multiple  frequency  and 
continuous  slip  vs.  stick-slip  or  stuck  behavior)  was  undertaken. 

It  was  demonstrated  earlier  that  the  coincidence  of  the  natural  frequencies  of  the 
primary  and  secondary  systems  is  a  particularly  critical  condition  since  it  leads  to  a 


maximum  transfer  of  energy  from  the  primary  to  the  secondary  where  it  can  be  dissipated 
by  friction.  The  effectiveness  of  the  frictional  device  is  thus  maximum  in  that 
configuration.  Accordingly,  the  present  study  of  the  effects  of  the  damping  coefficient  c, 
was  performed  under  this  matching  requirement,  i.e.  with  £,  =0.05.  Finally,  the 
coefficients  of  friction  were  selected  to  take  on  the  representative  values  jus  =juD  =0.25 
and  the  van  der  Pol  parameter  was  chosen  as  the  intermediate  value  of  8  =  0.025 .  The 
steady  state  response  of  the  system  of  Fig.  2.1  was  then  computed  for  a  series  of  values  of 
c,  and  c  and  the  largest  response  (e.g.  the  amplitude  of  the  motion  in  “single”  frequency 
solutions)  as  well  as  the  character  of  the  solutions  (i.e.  single  vs.  multiple  frequency, 
continuous  slip  vs.  stick-slip  vs.  stuck  configuration)  was  recorded,  see  Fig.  2.20  and 
2.21. 

It  should  first  be  noted  from  Fig.  2.20  the  existence  of  a  minimum  value  of  the 
largest  response  of  the  primary  mass  at  some  intermediate  value  of  the  damping 
coefficient  cv  When  c,  is  less  than  this  threshold,  it  would  appear  that  the  friction  and  the 
dashpot  reinforce  each  other  to  provide  an  increased  dissipation  and  correspondingly  a 
reduced  response.  A  more  refined  perspective  on  this  effect  may  be  obtained  from  the 
classification  of  the  various  solutions  involved,  see  Fig.  2.21.  Specifically,  it  is  seen  that 
there  exists  a  thin  layer  of  values  of  c{  near  c,=  0  in  which  the  character  of  the  solution 
changes  rapidly,  from  multiple  to  single  frequency.  In  fact,  a  comparison  of  Fig.  2.20  and 
2.21  indicates  that  it  is  also  in  this  small  zone  that  the  largest  response  decreases. 

It  is  thus  concluded  that  this  decrease  goes  hand  in  hand  with  a  change  of 
character  of  the  solution,  from  multiple  to  single  frequency.  To  complete  the  analysis,  the 
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time  histories  of  the  response  of  the  primary  system  (mass  M\ )  were  inspected  for  values 
of  cx  in  this  thin  layer,  i.e.  for  c,  e  [0,0.001],  see  Fig.  2.22.  It  is  clearly  seen  from  these 
figures  that  the  decrease  of  the  largest  response  is  in  fact  associated  with  a  decrease  of  the 
level  of  beating,  not  with  a  noticeable  decrease  in  the  mean  of  the  steady  state  response.  It 
remains  then  to  explain  how  this  change  of  character  takes  place  and  why  it  is  so 
sensitive  with  respect  to  c,. 

In  this  context,  consider  the  plots  of  the  largest  response  and  the  corresponding 
classification  obtained  with  c,  =0,  see  Fig.  2.14  and  2.15,  and  note  that  the  multiple 
frequency  solutions  only  exist  in  a  narrow  zone  of  stiffnesses  kx  around  kx  =0.05.  It  was 
argued  earlier  that  the  appearance  of  the  beating  is  associated  with  a  Hopf  bifurcation  in 
which  an  eigenvalue  of  the  perturbed  problem  reaches  the  unit  circle.  Damping  naturally 
provides  decay  over  a  given  time  and  thus  forces  these  eigenvalue  away  from  the  unit 
circle,  toward  the  origin.  For  c,  large  enough,  i.e.  0.001  in  the  case  considered,  the  largest 
eigenvalue  never  reaches  the  unit  circle  and  the  single  frequency  solution  is  maintained 
throughout.  Two  more  observations  can  be  made  in  this  regard.  First,  a  minimum  of  the 
largest  response  does  not  occur,  see  Fig.  2.14,  for  the  two  least  negative  values  of  c,  c  =  - 
0.02  and  -0.03,  consistently  with  Fig.  2.15  since  the  corresponding  solutions  for  c,  =0 
are  both  single  frequency.  Finally,  note  that  the  amplitude  of  the  single  frequency 
motions  for  cx=  0.001  and  kx-  0.05,  i.e.  6. 7-6.9  (see  Fig.  2.20),  are  very  close  to  the 
interpolation  of  the  values  obtained  from  Fig.  2.14  for  c,=  0  and  kx=  0.04  and  £,=  0.06 
suggesting  a  continuity  of  the  solution  as  a  function  of  in  the  absence  of  multiple 
frequency  bifurcation. 
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The  above  discussion  has  clarified  the  behavior  of  the  solutions  for  small  values 
of  q  but  has  not  addressed  the  increase  in  response  observed  in  Fig.  2.20  for  q  >  0.001. 
In  fact,  it  was  found  that  the  increase  in  damping  of  the  relative  motions  of  the  two  blocks 
leads  primarily  to  a  decrease  of  the  amplitude  of  these  motions  and  that  the  decrease  is 
large  enough  to  decrease  the  energy  dissipated  per  cycle,  even  with  the  increase  of  q . 
Thus,  the  amplitude  of  response  of  the  primary  mass  must  increase  as  noted  in  Fig.  2.20. 
The  decrease  of  amplitude  of  the  relative  motions  also  increases  the  likelihood  of  a 
sticking  and  thus  justifies  the  continuous  slip  to  stick-slip  transition  that  is  observed  on 
Fig.  2.21  for  all  c  values  as  q  increases.  For  very  large  values  of  q  ,  the  two  blocks  stick 
together  and  the  motion  is  governed  by  the  negative  dashpot  c  and  the  van  der  Pol 
restoring  force. 

It  was  desired  to  confirm  the  above  observations  in  another  situation,  namely  for  a 
different  value  of  the  van  der  Pol  coefficient  S .  To  this  end,  the  computations  were 
repeated  with  S  =  0  and  the  results,  largest  response  and  classification,  are  presented  in 
Fig.  2.23  and  2.24.  The  narrow  zone  in  which  both  transition  from  multiple  to  single 
frequency  and  decrease  of  the  largest  response  occur  is  again  present,  albeit  broader  in 
Fig.  2.23  and  2.24  as  compared  to  Fig.  2.20  and  2.21.  Further,  the  ensuing  increase  of  the 
amplitude  of  response  as  the  damping  coefficient  c,  is  increased  is  also  shown.  All 
results  thus  appear  consistent  with  those  of  Fig.  2.20  and  2.21. 
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Fig.  2.20  Largest  response  of  mass  Mx  as  a  function  of  the  damping  coefficient  c, 
for  different  values  of  the  damping  constant  c,  k .  =0.05,  <?=  0.025 ,  and  /j.  =  0.25 . 


SFCS  4  SFSS  ■  MFCS  *  MFSS 


Fig.  2.21  Domain  and  type  of  stable  solutions,  kx  =  0.05,  8 =  0.025,  //=  0.25 
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Fig.  2.25  Largest  response  of  mass  Mx  as  a  function  of  the  damping  coefficient  q 
for  different  values  of  the  damping  constant  c ,  A,  =  0.05 ,  5  -  0.025 ,  and  fx  =  0 . 

The  above  findings  and  discussion  demonstrate  the  benefit  of  a  small  viscous 
damping  in  addition  to  the  friction  as  a  mean  of  controlling  the  solution  character,  more 
notably  the  single  to  multiple  frequency  bifurcation.  A  larger  viscous  damping  is 
however  not  advisable  as  it  leads  to  a  reduction  of  the  dissipation  potential  of  the 
secondary  system. 

Having  established  the  role  in  friction  alone  and  in  concert  with  a  viscous 
damping  element,  it  is  desired  next  to  evaluate  the  stabilization  potential  of  the  viscous 
element  alone.  To  this  end,  the  computations  were  reiterated  but  for  p.  =  0  with  8  =  0.025 
and  the  corresponding  amplitudes  of  the  response  of  the  primary  mass  are  shown  in  Fig. 
2.25  as  function  of  q .  Interestingly,  it  is  found  that  the  steady  state  amplitude  is  zero  for 
a  finite  range  of  values  of  q,  q  e  [0.009, 0.018]  for  all  values  of  c  considered, 


45 


demonstrating  that  the  equilibrium  position  is  then  stable  even  with  the  presence  of  the 
negative  dashpot  c.  The  study  of  the  linearized  system  (i.e.  without  the  van  der  Pol  term) 
is  appropriate  to  explain  this  behavior  and  accordingly,  the  complex  mode  shapes  and 
eigenvalues  of  the  linear,  non-classically  damped  2-degree-of-freedom  system  were 
determined.  At  low  values  of  q ,  both  mode  shapes  involve  large  relative  motions  which 
further  increase  as  q  is  increased  so  that  both  modal  damping  ratios  increase  and 
eventually  become  positive  (e.g.  for  q  =  0.009  when  c  -  -0.12).  The  system  is  then  stable 
and  the  steady  state  amplitude  is  0.  Further  increases  in  cj  first  lead  to  small  increases  of 
the  relative  motions  of  both,  then  of  only  one  of  the  mode  shape  with  the  other  one 
exhibiting  a  decrease.  Notwithstanding  the  increase  in  q ,  this  reduction  in  relative 
motion  of  the  mode  shape  is  significant  enough  to  force  a  reduction  of  the  corresponding 
modal  damping  which  eventually  becomes  negative  (e.g.  forq  =0.018  when  c  =  -0.12). 
The  equilibrium  is  then  no  longer  stable.  Further  increases  of  q  lead  to  further 
reductions  of  the  relative  motions  of  that  mode  which  eventually  converges  to  [1  1]  as  q 
becomes  very  large.  That  is,  the  motions  associated  with  this  unstable  mode  correspond 
to  the  stuck  configuration  and  the  steady  state  response  of  the  system  is  governed  by  the 
van  der  Pol  restoring  force  on  the  stuck  system,  exactly  as  in  the  case  p  *  0. 
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Fig.  2.26  Largest  response  of  mass  Mx  as  a  function  of  the  damping  coefficient  c, 
for  different  values  of  the  damping  constant  c  and  coefficient  of  friction  // ,  noted  as 
c;ju .  kx-  0.05,  8  =  0.025.  The  crossings  between  the  corresponding  n  =  0  and 
H  =  0.25  curves  are  circled  in  red. 

A  comparison  of  the  curves  of  Fig.  2.20  and  2.25,  see  Fig.  2.26  for  a  few  values 
of  c,  is  of  particular  interest.  Most  importantly,  note  that  the  curves  corresponding  to 
H  =  0.25  and  //  =  0  cross  each  other  at  a  relatively  low  value  of  the  damping  coefficient 
ci .  For  values  of  q  smaller  than  this  threshold,  the  largest  response  of  the  steady  state 
motions  of  the  system  with  friction  is  less  than  its  frictionless  counterpart.  On  the 
contrary,  for  values  of  q  larger  than  the  threshold,  the  frictionless  system  is  the  one 
exhibiting  the  smallest  response.  In  fact,  in  the  interval  q  e  [0.009,0.018]  (for 
c>-0.12),  the  ffictionless  system  returns  to  equilibrium  while  the  one  with  friction 
maintains  a  fairly  large  amplitude  limit  cycle  oscillation!  In  this  case,  friction  may 
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significantly  increase  the  LCO  amplitude,  albeit  it  is  a  purely  dissipating  mechanism.  The 
explanation  of  this  surprising  observation  resides  in  the  sticking  potential  of  friction. 
Consider  for  example  the  response  of  the  system  with  cj  e  [0.009,0.018]  and  c>-0.12 
to  a  small  perturbation.  As  explained  above,  the  frictionless  system  has  a  stable 
equilibrium  for  these  parameter  values  and  thus  the  perturbation  decreases  with  time, 
returning  to  zero  after  an  infinite  time.  The  system  with  friction  does  not  exhibit  this 
behavior.  Specifically,  at  the  first  occurrence  of  a  zero  relative  velocity,  the  system 
becomes  stuck  because  the  forces  involved  (inertia,  dashpot  q  ,  and  spring  kj )  are  too 
small  to  maintain  slip.  Accordingly,  there  is  no  more  dissipation  and  the  response  of  the 
primary  mass  increases  rapidly  under  the  action  of  the  negative  dashpot.  Eventually,  the 
response  becomes  large  enough  to  force  slipping  and  dissipation  resumes  in  both  the 
dashpot  CJ  and  the  friction  element  and  convergence  to  a  limit  cycle  takes  place.  If  this 
limit  cycle  is  perturbed,  for  example  through  a  reduction  of  the  responses  of  the  two 
masses,  sticking  will  occur  at  the  next  occurrence  of  a  zero  relative  velocity  and  the 
motion  will  then  grow  back  under  the  sole  action  of  the  aerodynamic  effects  (negative 
dashpot). 

It  is  thus  concluded  that  friction  does  not  always  lead  to  a  decrease  of  the 
amplitude  of  LCO. 

2.4  Harmonic  Balance 

The  harmonic  balance  method  is  an  approximate  technique  for  the  estimation  of 
the  periodic  steady  state  response  of  nonlinear  systems.  In  this  approach,  the  response  is 
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assumed  in  the  form  of  a  limited  Fourier  series  the  unknown  constants  of  which  are 
evaluated  by  matching  as  many  cos  na>  t  and  sin  neat  components  as  possible  in  the 
equations  of  motion.  Note  in  the  present  context  of  self  excited  vibrations  that  the 
fundamental  frequency  co  is  also  unknown. 

The  harmonic  balance  method  will  be  used  here  to  obtain  an  approximation  of  the 
continuous  slip  response  of  the  system  of  Fig.  2.5.  At  the  contrary  of  the  exact 
formulation,  this  development  is  not  seriously  affected  by  the  value  of  k}  and  thus  was 
carried  out  directly  in  the  general  case  ■£  0  and  then  particularize  to  k,  =  0  for 
comparison  with  the  exact  formulation  and  the  numerical  simulation  results.  It  might 
seem  at  first  that  the  displacements  of  the  two  masses  should  both  be  expressed  as  sine 
and  cosine  of  all  frequencies  ncot  but  it  should  be  recognized  that  the  system  of  Fig.  2.5 
is  time  invariant  and  thus  that  any  specific  time  can  be  considered  as  the  temporal  origin 
t=  0.  It  is  particularly  convenient  here  to  select  this  point  to  match  a  maximum  value  of 
the  relative  displacement  of  the  two  masses.  Under  this  condition,  a  single  harmonic 
approximation  of  the  system  response  is 

Xj  =  Az  cos  co  t  +  Bz  sin  co  t  and  y  =  x 2  -  =  Y  cos  co  t .  (2-35) 

Next,  it  is  required  to  obtain  the  corresponding  representation  of  the  nonlinear 
forces,  i.e.  the  force  of  friction  here.  Given  the  above  form  of  the  relative  response,  it  is 
directly  found  that  a  one  term  Fourier  series  is 

/udN sgn(x)  «  4'--—  gin  (2-36) 


K 
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The  estimation  of  the  unknown  parameters  Az,  Bz,  Y,  and  co  is  then  achieved  by 
introducing  Eq.  (2-35)  and  (2-36)  in  the  equations  of  motion 

M{x\  +cx i  +  (k  +  k\ ) jcj  -k\X2  =  E31  and  M2X2  —k\X\  +&i*2  =--^31  (2-37) 

where 

F31  =//DArsgn(i2-^1).  (2-38) 

This  process  leads  to  a  set  of  four  nonlinear  algebraic  equations  for  the  unknowns. 
After  some  manipulations,  it  was  found  that 

(a-rp^-p^^-rp^C^P^rp^-rp2^  (2-39) 


B=- 


;  k  ’ 


=  {r+l)p2-14fiDN 
2  2  £rp3  Kk 


(2-40) 


and 


.  at.  .A 

kr-rp')  ' 


(2-41) 


where 

P  =  g>/g)'2;  r  =  M2/M{;  and  a=kl/k  (2-42) 

The  solution  of  Eq.  (2-39)  is  achieved  first  to  yield  the  value  of  the  frequency  p.  The 
remaining  parameters  are  then  estimated  from  Eqs  (2-40)  and  (2-41).  When  the  stiffness 
kx  vanishes,  the  cubic  characteristic  equation  (Eq.  2-39)  is  easily  solved  to  yield 
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A  two-frequency  harmonic  balance  approximation  was  also  obtained  by  considering  both 
cot  and  3 cot  sine  and  cosine  terms  in  both  and  y(t ).  That  is, 

jt]  =  A\  cos co t  +  B\  sin co t  +  A 3  cos 3 cot +  B^  sin 3 co t  (2-44) 


and 

y  -x 2  -  jcj  =  Cj  cos  co  t  +  D]  sin  co  t  +  C3  cos  3<a  t  +  D3  sin  3co  t .  (2-45) 

Since  the  temporal  origin  *=0  is  selected  to  correspond  to  a  peak  of  the  relative  response, 
one  must  have  y( 0)  =  0  or 


3  D3  +  D}  =  0 .  (2-46) 

Further,  the  2-harmonic  Fourier  series  of  the  friction  force  is  easily  shown  to  be 

/jdN  sgn(jr)  «  4 ——  sin  co  t  +  ———sin  3 co  t  .  (2-47) 

jc  In 

Then,  introducing  Eqs  (2-44),  (2-45),  and  (2-47)  in  the  equations  of  motion,  Eqs  (2-37) 
and  (2-38)  and  comparing  terms  in  sin  rut,  cosryi,  sin3cot  and  cos  3  cot  leads  to  the 
following  8  equations 


(l-jp2)^1+2C2p51-aC1=0 

(2-48) 

~2£iP  A  +(l  -  p2]^  -  a  D\  =  -  — 

(2-49) 

(l-9/72)^+6^53-crC3=0 

(2-50) 

■6apA3+{\-9p2)B3-<7D3=-~^ 

(2-51) 

-r  p1  A\  +  (cr  -  r  p2)C\  =0 

(2-52) 

-rp2  Bi  +(cr-rp2)Dl  =^- 

(2-53) 
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-9rp2A3+(<r-9rp2)c2=0  (2-54) 

-9rp2  B3+(a-9rp2)D2=^-  (2-55) 

Equations  (2-46)  and  (2-48)-(2-55)  represent  a  set  of  9  equations  for  the  unknown 
parameters  of  the  response  and  the  frequency  p.  The  solution  strategy  proceeds  as 
follows.  First,  Eqs  (2-49),  (2-51),  (2-53),  and  (2.55)  are  combined  to  produce  3 
homogenous  equations  and  1  non-homogenous  one.  The  former  equations  form  with  Eqs 
(2-46),  (2-48),  (2-50),  (2-52)  and  (2-54)  a  set  of  8  homogenous  equations  which  are  linear 
in  the  8  coefficients  A3 , ...,  D2  for  a  given  value  of p.  To  obtain  a  non-trivial  solution,  it 

is  thus  required  that  p  be  chosen  so  that  the  determinant  of  the  matrix  of  coefficients  of 
these  equations  vanish.  The  nonlinear  equation  resulting  from  this  condition  is  a 
polynomial  of  order  16  mp.  For  each  of  its  roots,  the  coefficients  A\ , ....  D2  can  then  be 

evaluated  by  solving  the  system  of  linear  equations  formed  by  7  of  the  8  homogenous 
equations  and  the  non-homogenous  one. 

The  above  procedures  were  followed  to  obtain  1-  and  2-harmonic  approximation 
of  the  continuous  slip  response  of  the  system  of  Fig.  2.5  and  their  accuracy  was  assessed 
by  comparison  with  the  numerical  simulation  and/or  the  exact  formulation.  The  special 
case  =  0  was  considered  first  and  the  results  are  presented  in  Fig.  2.4.  Clearly,  there 
are  two  frequencies  satisfying  Eq.  (2-43)  (one  for  the  +  sign,  the  other  for  the  -  sign)  and 
thus  the  harmonic  balance  solution  yields  two  estimates  of  the  continuous  slip  solution, 
one  of  which  (R)  approximates  the  true  response  while  the  other  is  a  fictitious  solution 
(F).  It  is  seen  from  Fig.  2.4  that  the  accuracy  of  the  harmonic  balance  (solution  R)  is  not 
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very  good  in  the  small  domain  of  existence  of  these  continuous  slip  solutions.  Further,  the 
inclusion  of  a  second  harmonic  does  not  seem  to  improve  substantially  the  accuracy  of 
the  estimates  of  the  amplitude  and  frequency  of  response,  see  Fig.  2.4  in  which  (HBM)- 
CS  and  (HBM2)-CS  are  the  1-  and  2-harmonic  approximations  of  the  continuous  slip 
solutions. 


Fig.  2.27  Relative  errors  (in  %)  in  the  maximum  response  of  mass  Mx  as  predicted 
by  the  harmonic  balance  method  with  1  and  2  harmonics  as  functions  of  kx , 

M2  =  0.05,  c  =  -0.07  (Note  the  split  horizontal  axis). 

The  accuracy  of  the  harmonic  balance  method  was  also  assessed  for  a  broad  range 
of  ki  *  0  situations.  Shown  in  Fig.  2.27  are  typical  relative  errors  on  the  maximum 
response  obtained  in  the  SFCS  domain  by  using  the  1-  and  2-harmonic  approximations 
and  corresponding  to  a  frequency  close  to  the  one  obtained  by  numerical  simulation. 
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Given  these  good  qualitative  results,  the  harmonic  balance  method  was  also  used  to 
estimate  the  boundary  of  the  stabilization  domain,  see  the  curves  HBM1F  and  HBM2F  on 
Fig.  2.6.  It  appears  from  all  these  results  that  the  harmonic  balance  estimates  were 
qualitatively  correct  in  predicting  the  amplitude  of  the  response  and  that  the  addition  of  a 
second  harmonic  provided  an  improved  quantitative  matching  with  the  simulation  results. 

The  difference  in  accuracy  of  the  harmonic  balance  method  for  ftj  =  0  and  k\  *  0 
may  result  from  the  naturally  harmonic  character  of  the  response  of  the  sliding  mass  in 
the  latter  case  as  opposed  to  the  parabolic  time  histories  obtained  in  the  former  case. 


CHAPTER  3 


DESIGN  AND  TESTING  OF  A  FRICTION  DEVICE 

The  second  part  of  this  thesis  focuses  on  the  design,  fabrication,  and  testing  of  a 
friction  device  similar  to  the  disk-torsional  spring  system  of  Fig.  1.2  which  was  later 
tested  in  the  Deutsches  Zentrum  fur  Luft  und  Raumfahrt  (DLR)  transonic  wind  tunnel. 
The  results  of  Chapter  2  serve  here  as  stepping  stone  for  the  parameter  selection  of  the 
device  components. 

3.1  Design  and  Fabrication 

The  initial  plan  for  the  friction  device  was  an  installation  inside  the  airfoil, 
consistently  with  Fig.  1.2,  but  this  option  led  to  some  daunting  challenges,  e.g. 
fabricating  a  novel  airfoil  (cost  constraint),  fabricating  a  novel  airfoil  with  enough  space 
internally  for  the  device,  having  access  to  the  device  once  the  airfoil  was  finished,  etc.  In 
this  light,  it  was  proposed  by  DLR  to  mount  the  device  rigidly  to  the  airfoil  but  outside  of 
the  wind  tunnel,  i.e.  on  the  balance  rigs  located  on  either  side  of  the  wind  tunnel  and 
airfoil.  The  first  bending  natural  frequency  of  the  airfoil  is  known  to  exceed  100Hz  while 
the  vibrations  of  the  device  were  expected  to  take  place  near  the  flutter  frequency,  i.e. 
around  30Hz.  It  was  thus  argued  that  the  airfoil  would  appear  rigid  and  thus  that  the  2- 
dimensional  character  of  the  problem  would  not  be  compromised  by  the  friction 
mechanism  being  lumped  at  the  ends  of  the  airfoil  as  opposed  to  being  distributed  along 
the  span. 
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Even  then,  the  design  of  the  friction  device  was  very  substantially  constrained  in 
regards  to  weight,  size,  and  to  a  smaller  extent  cost.  From  the  standpoint  of  validation  of 
the  present  findings,  it  was  desired  to  achieve  an  inertia  ratio  of  at  least  5%  for  the  effects 
of  friction  to  be  most  noticeable  but  to  maintain  the  occurrence  of  the  flutter/LCO  in  the 
transonic  range,  it  was  necessary  to  limit  this  ratio  to  about  3%.  The  size  limitations 
originated  from  the  need  to  maintain  the  balance  rig  in  its  present  configuration. 

The  lack  of  access  to  the  balance  rig,  and  thus  to  the  friction  devices,  during  the 
wind  tunnel  testing  was  the  next  important  issue  to  be  resolved.  Indeed,  to  minimize  the 
forces  on  the  wind  tunnel  walls  and  to  permit  the  use  of  porous  walls,  the  central  section 
of  the  wind  tunnel,  including  the  balance  rigs,  is  confined  to  a  plenum  which  is 
depressurized,  and  thus  not  accessible,  during  the  operation  of  the  wind  tunnel.  Further, 
repressurization  is  possible  when  the  wind  tunnel  is  down  but  the  repressurization/ 
depressurization  process  requires  about  1  hour,  i.e.  a  significant  portion  of  the  total 
testing  time  available.  For  validation  of  the  theoretical  results,  it  was  however  desired  to 
proceed  with  a  parametric  study  of  the  effects  of  the  disk  properties  (inertia  I q  and/  or 
stiffness  kD)  and  friction  force  (juN)  on  the  response  of  the  airfoil-device  system.  This 
situation  motivated  the  search  for  ways  to  vary  the  inertia/stiffness  of  the  disk  and  the 
friction  force  remotely.  While  no  simple  means  for  changing  the  inertia/  stiffness 
remotely  was  devised,  it  was  decided  to  vary  the  friction  force  by  regulating  the  normal 
force  through  the  compression  of  linear  springs  by  a  linear  actuator. 

While  the  use  of  linear  actuators  provided  the  remote  variation  capability  desired, 
it  also  incurred  a  large  increase  in  dead  weight  incompatible  with  the  tight  limits  of  mass 
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and  inertia  to  maintain  the  LCO  in  the  transonic  range.  The  only  solution  found  to  this 
problem  was  the  mounting  of  the  actuator  system  on  the  laboratory  fixed  frame  of  the 
balance  rig.  The  solution  of  the  weight  issue  then  triggered  a  new  problem:  the 
transmission  of  the  normal  force  and  normal  force  only  (no  tangential  force  nor  moment) 
from  the  laboratory  fixed  actuator  system  to  the  disk  moving  in  the  plane  perpendicular  to 
the  force  to  be  applied.  This  last  significant  difficulty  was  resolved  primarily  through  the 
use  of  the  cage  of  a  thrust  ball  bearing  inserted  between  the  disk  and  the  actuator  system. 
The  well  lubricated  balls  of  the  bearing  allowed  the  relative  motion,  transferred  the 
normal  force,  and  generated  only  a  minimal  side  (tangential)  force  with  good  lubrication. 
Finally,  the  offset  of  the  normal  forces  on  the  actuator  and  the  disk  induced  by  the  motion 
of  the  latter  would  normally  generate  a  moment  on  the  actuator  system  (severely  limited 
by  the  manufacturer)  and  potentially  (through  any  deformation)  a  non-uniformity  of  the 
pressure  applied  on  the  disk.  Both  of  the  potential  issues  were  resolved  by  inserting  a 
spherical  ball  contact  along  the  path  of  the  force  transmission  in  the  actuator  system.  That 
is,  the  actuator  system  pushed  on  the  ball  which  then  in  turn  pushed  on  the  thrust  bearing 
cage.  The  point  contacts  between  the  ball  and  the  contacting  shaft  insured  the 
transmission  of  a  normal  force  only. 

The  detailed  design  of  the  friction  device  and  actuator  system  was  then 
undertaken,  with  regular  interaction  with  DLR,  within  the  software  SolidWorks  using  an 
accurate  CAD/CAM  model  of  the  balance  rig  provided  by  DLR.  The  final  detailed 
model,  shown  in  Fig.  3. 1-3.3,  can  be  summarized  as  follows.  The  torsional  frictional 
device  is  composed  of  a  disk  connected  to  the  airfoil  (outside  of  the  wind  tunnel  section) 
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through  a  torsional  spring  machined  to  a  specific  stiffness.  The  disk  is  sandwiched 
between  two  steel  rings  providing  the  friction.  The  normal  force  is  provided  by  the 
actuator  system  shown  on  the  right  in  Fig.  3.3,  i.e.  generated  by  a  linear  actuator  based  on 
a  step  motor.  The  motor  pushes  on  three  parallels  springs  to  generate  a  controllable  force 
the  magnitude  of  which  is  monitored  through  a  load  cell  mounted  in  series  and  a  LVDT 
recording  the  translation  of  the  motor.  The  force  is  then  transmitted  from  the  load  cell  to 
a  small  spherical  ball,  thereby  preventing  the  transmission  of  anything  else  than  a  normal 
force,  and  then  to  a  “thrust  washer”.  This  subsystem  is  composed  of  a  disk  and  the  cage 
of  a  thrust  bearing,  the  latter  of  which  is  squeezed  between  the  thrust  washer  disk  and  the 
external  friction  ring  of  the  device.  The  bearing  cage  permits  the  transmission  of  the 
normal  force  even  with  relative  motion  of  the  two  disks  (thrust  washer  and  friction  disks). 
Note  that  the  bearing  cage  was  loosely  held  in  place  (preventing  it  from  falling)  by  three 
springs  connected  to  long  arms  themselves  attached  to  the  support  of  the  actuator  system. 
The  actuator  permits  the  variation  of  the  normal  force  while  a  change  of  the  natural 
frequency  of  the  frictional  device  is  allowed  by  the  addition/removal  of  masses  at  the 
periphery  of  the  disk  and  through  the  switch  of  torsional  springs  (2  different  stiffnesses 
were  available). 

The  fabrication  of  the  various  parts  was  performed  by  the  two  ASU  machine 
shops  as  well  as  an  outside  machine  shop. 


Fig.  3.1  Exploded  view  of  the  disk,  spring,  and  friction  ring  system. 


Fig.  3.2  Complete  view  of  the  airfoil,  friction  device,  and  actuator  system  in  the 
balance  rig  in  the  DLR  wind  tunnel. 


Fig.  3.3  Complete  view  of  the  friction  device  and  actuator  system  in  ASU 

vibrations  laboratory. 

3.2  Instrumentation 

The  instrumentation  of  the  DLR  wind  tunnel  permits  the  measurement  of  the 
pressure  at  a  series  of  locations  on  the  wind  tunnel  walls,  of  the  lift  and  moment  on  the 
airfoil,  of  the  plunging  displacement  and  pitching  angle,  and  of  the  displacements  at 
selected  locations  on  the  airfoil.  Of  primary  importance  for  the  present  validation  were 
the  plunging  displacement  and  pitching  angle  but  also  the  relative  angle  of  torsion  9.  This 
measurement  was  obtained  through  the  reading  of  the  strain  in  the  coils  of  the  torsional 
spring  (from  a  strain  gauge  mounted  directly  on  the  thin  coils).  The  calibration  factor 
between  angle  and  strain  was  obtained  experimentally  for  each  spring.  The  procedure  for 
exciting  and  controlling  the  airfoil  and  other  details  on  the  experimental  capabilities  were 
the  same  as  those  used  in  prior  investigation  of  transonic  LCO  at  DLR. 

3.3  Testing  in  ASU  Vibrations  Laboratory 

The  primary  purpose  of  the  testing  that  took  place  in  the  ASU  Vibrations 
Laboratory  was  to  validate  the  design  and  fabrication,  i.e.  to  demonstrate  that  relative 
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motions  of  the  disk  with  respect  to  its  support  could  take  place  and  that  these  motions 
could  be  controlled  by  the  force  actuator.  In  particular,  it  was  desired  to  assess  the 
capability  of  the  device  to  exhibit  stick,  stick-slip,  and  continuous  slip  motions.  To  this 
end,  the  device  was  mounted  on  a  bar  of  a  design  similar  to  the  device  support  in  the 
wind  tunnel.  This  bar  was  then  attached  to  a  frame,  itself  mounted  on  the  Ling  shaker  of 
the  ASU  Vibrations  Laboratory  as  to  impart  plunging  and  mostly  pitching  of  the  bar.  The 
shaker  was  used  to  provide  a  harmonic  excitation  at  or  near  30Hz,  i.e.  similar  to  what 
could  be  expected  in  the  wind  tunnel.  The  normal  force  applied  was  varied  to  assess  the 
capability  of  the  device  to  exhibit  the  required  motions. 

The  exact  classification  of  the  observed  motions  was  not  always  very  easy 
because  of  the  noise  present  in  the  measurement,  about  10  microstrains.  Nevertheless,  it 
is  believed  that  all  three  types  of  motions  were  indeed  noted.  Stuck  motions  were 
observed  for  low  excitation  levels  from  the  shaker  and  large  normal  forces  as  expected 
and  are  not  discussed  further  here.  Continuous  slip  responses  were  initially  noted  at  zero 
normal  force,  see  Fig.  3.4,  but  also  for  nonzero  forces,  see  Figs  3.5  and  3.6.  Stick-slip 
motions  were  observed  first,  see  Figs  3.7  and  3.8,  with  sticking  taking  place  at  the  peak 
displacements  but  an  intermediate,  asymmetrical  sticking  phase  has  also  been 
encountered,  see  Figs  3.9  and  3.10.  It  is  unclear  what  created  the  asymmetry  in  the 
response,  two  possible  sources  are  the  torsional-axial  coupling  of  the  spring  and/or  the 
variations  of  the  normal  force  (see  the  load  cell  time  history  on  Figs  3.9  and  3.10)  which 
could  have  resulted  from  vibrations  of  the  shaker.  In  regards  to  the  former  issue,  note  that 
the  spring  wants  to  displace  axially  as  it  is  twisted  in  a  direction  (compression  or 


extension)  that  is  dependent  on  the  direction  of  torsion  (clockwise  or  counterclockwise). 
If  the  force  applied  by  the  friction  rings  is  not  the  same  on  the  two  sides,  this  axial-torsion 
coupling  would  induce  a  variation  of  the  normal  force  that  is  related  to  the  torsion  angle 
which  might  produce  the  asymmetrical  behavior  observed. 


Fig.  3.4  Continuous  slip  response  with  zero  normal  force  at  30Hz  excitation. 


Fig.  3.6  Continuous  slip  response  with  nonzero  normal  force  at  30Hz  excitation. 


Fig.  3.8  Stick-slip  response  with  nonzero  normal  force  at  30Hz  excitation. 


Fig.  3.10  Stick-slip  response  with  nonzero  normal  force  at  30Hz  excitation. 
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3.4  Testing  in  the  DLR  Transonic  Wind  Tunnel 

Albeit  complex,  the  friction  device/actuator  system  combination  worked  perfectly 
as  designed  and  led,  for  the  first  time,  to  the  generation  and  recording  in  the  DLR  wind 
tunnel  of  transonic  limit  cycles  with  active  internal  friction  at  Mach  numbers  close  to  0.6. 

The  first  important  observation  to  be  drawn  is  that  such  experimental  analyses  are 
possible  and  that  internal  relative  motions  can  indeed  be  excited  in  LCO  scenarios. 

Some  of  the  recorded  LCO  are  shown  in  Fig.  3.11  and  3.12.  Although  only  a 
limited  set  of  LCO  records  were  obtained,  some  observations  can  nevertheless  be  made. 
Consider  for  example  the  evolution  of  the  time  histories  of  the  response  (plunging, 
pitching,  and  torsion)  as  a  function  of  the  Mach  number.  While  the  plunging  and  pitching 
time  histories  show  a  clean  single  frequency  behavior  seemingly  independent  of  the 
Mach  number,  the  relative  torsion  angle  displays  a  much  richer  behavior.  Indeed,  for 
Ma  =  0.598  and  neglecting  the  high  frequency  noise,  the  time  history  of  9  appears  at 
first  to  be  almost  purely  a  single  frequency  solution.  A  closer  inspection  however 
suggests  that  the  rise  of  9  is  steeper  than  its  fall,  especially  near  the  peak.  As  the  Mach 
number  is  increased,  this  effect  becomes  more  pronounced  and  a  flattened  peak  is 
obtained  for  Mv> 0.606,  suggesting  the  possibility  of  a  sticking  phase.  The  lack  of 
symmetry  of  this  flattening  implies  a  full  Fourier  series,  including  even  harmonics  at  the 
contrary  of  the  solutions  obtained  computationally  so  far. 
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Fig.  3.11,  continued 
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e)  Mx  =  0.606  f)  Mm  =  0.608 

Fig.  3.12  Zoomed  time  histories  of  the  plunging  displacement  ("heave"),  the  pitching 
angle  ("alpha"),  and  the  relative  torsional  angle  ("alWD")  recorded  in  the  DLR  wind 

tunnel  for  a  normal  force  of  21  N. 
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a)  Overall  time  histories  at  a  normal  force  of  50  N 


b)  Zoomed  time  histories  at  a  normal  force  of  50  N 

Fig.  3.13  Time  histories  of  the  plunging  displacement  ("heave"),  the  pitching  angle 
("alpha"),  and  the  relative  torsional  angle  ("alWD")  recorded  in  the  DLR  wind 

tunnel  for  =  0.608 . 
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c)  Repeatability  check  at  50  N 
Fig.  3.13,  continued 

What  about  the  effects  of  the  normal  force?  Shown  in  Fig.  3.13(a)  and  (b)  are  the 
time  histories  of  the  three  variables  for  a  Mach  number  of  0.608  and  a  normal  force  of 
50N,  as  opposed  to  the  21N  used  in  Fig.  3.11(f)  and  3.12(f).  A  comparison  of  the  time 
histories  of  the  2  IN  and  50N  is  given  in  Fig.  3.13(b)  from  which  it  is  concluded  that  the 
increase  in  normal  force  has  led  to  a  decrease  of  the  torsional  response  and, 
correspondingly,  to  a  small  increase  in  the  plunging  and  pitching  amplitudes.  The  small 
changes  that  are  involved  in  Fig.  3.13(b)  ought  to  be  confirmed  as  genuine  and  not  the 
result  of  noise  or  variability.  To  this  end,  several  cases  were  repeated  to  provide  a 
baseline  for  the  assessment  of  the  noise  and  variability.  In  fact,  a  repeat  measurement  of 
the  LCO  for  A/  =  0.608  and  a  50N  normal  force  is  available  and  the  two  sets  of  time 

«3 

histories  are  plotted  in  Fig.  3.13(c).  It  is  clear  from  this  figure  that  the  repeatability  is 
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high  so  that  small  differences,  such  as  those  noted  in  Fig.  3.13(b)  can  indeed  be 
interpreted  as  genuine  changes  in  the  response. 

In  closing  this  chapter,  it  should  be  noted  that  the  flutter  speed  of  the  system  in 
stuck  mode  was  found  to  be  higher  than  =  0.608  in  an  earlier  set  of  tests.  In  fact,  it 
was  necessary  to  excite  the  airfoil  by  a  shaker  during  the  tests  to  obtain  the  LCO 
measurements  of  Figs  3.1 1-3.13.  It  is  thus  concluded  from  these  two  observations  that  the 
results  presented  are  of  a  subcritical  nature. 


CHAPTER  4 


CONCLUSIONS 

The  focus  of  this  thesis  has  been  on  the  assessment  of  the  effects  of  internal 
friction  on  the  response  of  aeroelastic  systems  exhibiting  either  explosive  flutter  or  limit 
cycle  oscillations.  In  the  first  part  of  this  work,  a  two-degree-of-freedom  system  has  been 
extensively  studied  that  models  a  wing  with  an  internal,  sliding  component.  The  largest 
mass  of  the  system  represents  the  wing  while  the  smallest  one  is  associated  with  the 
sliding  component.  A  dashpot  of  negative  constant  was  introduced  to  model  the  unstable 
linear  aerodynamic  effects,  while  an  additional  van  der  Pol  restoring  force  was  included 
to  characterize  possible  stable  nonlinear  aerodynamic  effects.  Both  of  these  components 
were  assumed  to  act  only  on  the  largest  mass  (wing).  The  sliding  mass  was  finally 
assumed  to  be  connected  to  the  wing  through  a  friction  slider,  a  spring  ( k{  >:  0 ),  and 
possibly  a  dashpot  (positive  viscous  damping,  c\  >0). 

The  system  was  analyzed  first  without  a  van  der  Pol  term.  Further,  in  the  absence 
of  elastic  deformation  during  sliding,  i.e.  for  =  0 ,  it  was  shown  that  an  exact 
procedure  can  be  followed  to  transform  the  nonlinear  equations  of  motion  into  a  set  of 
nonlinear  algebraic  equations  and  that  a  stabilization  of  the  unstable  aerodynamic  forces 
was  possible.  The  most  negative  damping  that  could  be  stabilized  in  this  configuration 

was  in  fact  shown  to  be  2  r/n  (for  r  «  1)  where  r  denotes  the  ratio  of  the  two  masses 
of  the  system.  The  system  was  found  to  exhibit  periodic  motions  with  possible  stick 
phases.  The  introduction  of  a  nonzero  spring  constant  was  shown  to  lead  to  an  increased 
stabilization  potential  that  is  maximum  in  a  “tuned  damper”  configuration,  i.e.  when  the 
natural  frequency  of  the  primary  component  (the  wing  alone)  closely  matches  the  natural 
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frequency  of  the  secondary  alone  (sliding  mass  on  spring  k\ ).  In  addition  to  this  benefit, 
a  broader  set  of  motions  was  also  observed.  Besides  periodic  (referred  to  as  single 
frequency  here)  solutions,  multiple  frequency  motions  were  also  found  to  exist, 
especially  near  the  tuned  damper  configuration  while  the  single  frequency  solutions  were 
typically  obtained  near  the  stability  border.  Both  continuous  slip  and  stick  slip  multiple 
frequency  motions  were  noted  which  were  shown  to  be  either  aperiodic  or  chaotic. 
Further,  the  transition  from  single  to  multiple  frequency  was  shown  to  be  a  Hopf 
bifurcation.  It  was  finally  observed  that  the  maximum  steady  state  response  of  the  wing  is 
typically  minimum  for  the  single  frequency  motions  occurring  just  next  to  the 
bifurcation.  A  sharp  increase  in  response  level  is  obtained  after  bifurcation  that  results 
from  the  beating  induced  by  the  presence  of  multiple  frequencies. 

The  inclusion  of  a  van  der  Pol  term  did  not  significantly  change  the  above 
findings  with  the  noted  exception  that  the  system  is  always  stable  even  in  the  absence  of 
friction  in  which  case  it  experiences  a  continuous  slip  single  frequency  motion.  In  fact,  it 
was  observed  that  the  van  der  Pol  term  leads  to  a  significant  reduction  of  the  large 
amplitude  beating  excursions  obtained  otherwise  in  connection  with  multiple  frequency 
solutions  and  appears  to  postpone  the  Hopf  bifurcation.  In  regards  to  stick  slip  vs. 
continuous  slip,  the  decrease  in  response  implied  by  the  presence  of  the  van  der  Pol  term 
was  found  to  increase  the  likelihood  of  stick  slip  motions  which  typically  occur  at  lower 
response  levels  than  continuous  slip  solutions  do.  The  effects  of  the  coefficients  of 
friction  (the  static  and  dynamic  coefficients  were  assumed  equal  here)  was  studied  next. 
It  was  first  noted  that  the  amplitude  of  continuous  slip  motions  (single  or  multiple 
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frequency)  is  linearly  scaled  by  the  coefficient  of  friction  and  this  parameter  plays  no 
further  role  in  these  cases.  When  stick  phases  occur,  a  more  complex  dependency  on  the 
coefficient  of  friction  is  obtained.  Next,  it  was  found  that  2  stable  solutions  can  be 
obtained  for  low  coefficients  of  friction.  One  such  solution  has  a  low  amplitude  and  is 
dominated  by  friction  effects  alone  with  little  influence  of  the  van  der  Pol  term.  The 
reverse  statement  holds  for  the  other  solution,  the  amplitude  is  large  with  the  stabilization 
dictated  by  the  van  der  Pol  force  with  little  effect  of  friction.  As  the  coefficient  of  friction 
increases,  these  two  solutions  become  closer  together  and  eventually  merge.  The 
amplitude  of  response  then  appears  to  increase  monotonically  with  further  increases  of 
the  coefficient  of  friction.  The  effects  of  a  second,  viscous  damper  in  parallel  to  the 
friction  element  was  finally  considered.  It  was  found  that  the  two  damping  mechanisms 
do  not  necessarily  reinforce  each  other  and  that  the  friction  hinders  the  viscous 
dissipation  when  stick  occurs.  In  fact,  nonzero  amplitude  limit  cycle  oscillations  were 
noted  with  friction  in  situations  where  the  equilibrium  (zero  response)  would  in  fact  be 
stable  without  friction. 

The  second  part  of  the  thesis  focused  on  the  design,  fabrication,  and  testing  of  an 
actual  friction  device  on  an  airfoil.  The  basic  design  of  the  system  is  a  rotating  disk 
connected  to  the  airfoil  by  a  torsional  spring  of  natural  frequency  closely  matching  that  of 
the  airfoil  to  achieve  the  tuned  damper  arrangement  discussed  above.  The  spring  was 
instrumented  with  strain  gauges  to  measure  the  angle  of  torsion  of  the  system.  Finally,  a 
complex  system  was  designed  to  produce  the  required  normal  force  to  induce  friction  on 
the  moving  airfoil  without  transmitting  any  shear  or  moment.  The  preliminary  testing  of 
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the  system  was  accomplished  in  shaker  tests  conducted  in  the  ASU  Vibrations 
Laboratory.  The  system  was  shown  to  exhibit  the  basic  motions  expected,  i.e.  continuous 
slip  and  stick  slip  responses.  The  friction  device  was  then  installed  on  a  NLR  7301  airfoil 
in  the  DLR  transonic  wind  tunnel  operating  at  a  Mach  number  near  0.6.  Subcritical  limit 
cycle  oscillations  were  observed  in  which  the  friction  played  a  significant  role.  The 
limited  data  obtained  demonstrates  a  slight  decrease  in  response  by  increasing  the  normal 
force  applied,  i.e.  by  increasing  the  level  of  friction. 
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ABSTRACT 


The  appearance  of  limit  cycle  oscillations  in  aeroelastic  systems  has  usually 
been  associated  with  nonlinearity  in  the  aerodynamics  and/or  in  the  structural  restoring 
forces.  It  has,  however,  been  recently  suggested  that  nonlinearity  in  the  damping 
mechanism,  more  notably  friction  between  a  small  moving  part  (or  many  such  parts  such 
as  stores  and  missiles)  and  the  wing,  may  indeed  be  a  source  of  post-flutter  limit  cycle 
oscillations.  The  present  work  provides  a  numerical  validation  of  this  expectation  by 
studying  the  response  at  and  above  flutter  of  a  flat  plate  and  a  NACA0012  airfoil  both 
exhibiting  an  internal  friction  mechanism  and  placed  in  a  inviscid  and  incompressible 
flow. 

The  aerodynamic  of  the  flat  plate,  dictated  by  the  Theodorsen  function,  is  linear 
and  thus  the  frictionless  system  is  either  stable  or  flutters.  Time  marching  computations 
were  carried  out  with  a  rational  approximation  of  the  Theodorsen  function  on  the  system 
with  an  internal  friction  component.  These  results  demonstrate  the  existence  of  a  rich 
collection  of  limit  cycle  oscillations  in  a  range  of  flow  speeds  extending  above  the  flutter 
speed  of  the  flat  plate  with  its  frictional  device  stuck  by  about  3%  for  a  device  whose 
inertia  is  5%  of  the  inertia  of  the  flat  plate.  This  significant  increase  in  stability  zone 
demonstrates  the  good  potential  stabilization  of  friction. 

The  second  airfoil,  a  NACA0012  section  also  placed  in  an  inviscid, 
incompressible  flow,  had  been  found  in  earlier  studies  to  exhibit  limit  cycle  oscillations  at 
and  slightly  above  its  flutter  speed.  The  inclusion  of  a  friction  device  is  once  again  seen 
to  be  beneficial  as  significant  reductions  of  the  amplitude  of  the  limit  cycles  (e.g.  by  a 
factor  of  3)  can  be  obtained  by  appropriately  “tuning”  the  friction  device. 

iii 


The  nature,  i.e.  periodic,  aperiodic,  chaotic  and  continuous  vs.  stick  slip,  and 
magnitude  of  all  limit  cycle  oscillations  encountered  are  discussed  in  details. 
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CHAPTER  1 
INTRODUCTION 

Limit  Cycle  Oscillations  (LCO)  have  been  a  prevalent  aeroelastic  problem  on 
several  current  fighter  aircraft.  This  phenomenon  usually  occurs  for  aircraft  with  external 
stores  throughout,  but  not  limited  to,  the  transonic  flight  regime1,2,  although  a  business  jet 
wing  LCO  was  also  reported  recently3.  Complicated  by  the  problem  geometry,  e.g.  the 
aircraft-store  system,  the  LCO  mechanisms  still  remain  to  be  fully  understood.  In  fact, 
there  exist  few  analytical  techniques  available  for  LCO  prediction  and  an  insufficient 
understanding  of  its  physics. 

With  aerodynamic  feedback,  LCO  are  sustained  periodic  oscillations  which 
neither  increase  nor  decrease  in  amplitude  over  time  for  a  given  flight  condition.  A  series 
of  researchers,  notably  Cunningham  and  Meijer4,  believe  that  the  wing/store  LCO  is  a 
purely  aerodynamics  phenomenon,  largely  due  to  transonic  shock  oscillation  and  shock 
induced  flow  separation.  This  LCO  scenario,  which  is  referred  to  as  the  Transonic 
Shock/Separation  (TSS)  model,  has  been  suggested  by  Edwards  to  be,  with  viscous 
effects,  one  of  the  major  factors  contributing  to  transonic  LCO  for  wings3. 

In  1998,  Chen,  Sarhaddi  and  Liu2  offered  a  radically  different  LCO  model  based 
on  the  observation  that  wing/store  LCO  can  be  a  post-flutter  phenomenon  whenever  the 
flutter  mode  contains  low  unstable  damping.  This  type  of  flutter  mode  is  called  a  “hump 
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mode”.  Since  the  aircraft  structure  usually  exhibits  some  structural  nonlinearity,  such  as 
friction,  free-play,  etc.,  the  moderate  growth  of  amplitudes  corresponding  to  the  low 
negative  aerodynamic  damping  may  be  suppressed.  The  result  is  then  a  steady  state 
oscillation.  The  consideration  of  friction  in  wing/wing-store  LCO  studies  is  rather 
recent5.  Here  it  is  referred  to  as  the  nonlinear  structural  damping  (NSD)  model  of  the 
wing/store  LCO. 

It  should  be  noted  that  the  consideration  of  friction  in  flutter  analyses  is  not  new, 
as  it  is  recognized  in  the  turbomachinery  community  that  the  friction  at  the  blade  root  and 
in  shrouds  plays  a  definite  stabilizing  role.  In  fact,  this  observation  has  led  to  the  use  of 
specially  designed  friction-based  damping  systems,  referred  to  as  “underplatform” 
dampers,  to  damp  blade  vibrations  and  increase  the  flutter  speed  (e.g.  see  Ref.  6  and  7). 
Notwithstanding  these  similarities,  there  are  substantial  differences  in  the 
aerodynamic/structure  interactions  in  wing/wing-store  LCO  and  in  the  turbomachinery 
blade  flutter  that  render  each  problem  specific. 

In  light  of  the  above  discussion  the  focus  of  the  present  investigation  and  of  a 
companion  one,  Ref.  8,  is  on  providing  a  validation  of  the  potential  of  friction  to  stabilize 
an  unstable  aerodynamics  and/or  affect  the  magnitude  of  existing  limit  cycle  oscillation. 

The  purely  structural  dynamic  set  up  of  Fig.  1 . 1  was  a  adopted  in  Ref.  8  as  a 
simple  model  of  an  aeroelastic  system  near/above  flutter;  the  damping  coefficient  c  was 
taken  negative  to  simulate  small  amplitude  unstable  aerodynamic  effects  and  a  van  der 
Pol  restoring  force  could  be  included  to  represent  a  potential  stabilization  of  aerodynamic 
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origin.  In  this  fashion,  the  effects  of  friction  could  be  assessed  on  purely  fluttering 
systems  (without  the  van  der  Pol  term)  and  those  exhibiting  aero-driven  limit  cycles  (with 
van  der  Pol  term).  In  fact,  a  comprehensive  quantitative  analysis  of  this  system  was 
achieved,  most  notably  the  stabilization  potential  of  friction  was  substantiated,  and  the 
effects  of  the  various  parameters  clearly  described. 


Fig.  1.1  2-degree-of-freedom  system  with  van  der  Pol  (v.d.p)  force 

The  second  part  of  Ref.  8  described  the  design  of  an  actual  friction  device  and 
reported  the  results  of  wind  tunnel  tests  conducted  in  the  DLR  transonic  wind  tunnel. 
Some  of  the  finding  of  this  experimental  effort  confirmed  the  analysis  conducted  on  the 
system  of  Fig.  1.1  while  others  were  unexpected. 

In  this  light,  the  specific  goal  of  the  present  investigation  is  to  parallel  the  analysis 
of  Ref.  8,  but  on  actual  aeroelastic  systems  as  opposed  to  the  structural  dynamic  model  of 
Fig  1.1,  and  thus  to  bridge  the  two  components  of  this  earlier  study. 

In  selecting  an  appropriate  structural  model  for  a  full  aeroelastic  analysis,  it  is  first 
necessary  to  describe  in  details  how  friction  will  be  introduced.  The  representation  of  the 
friction  effects  is  particularly  challenging  as  it  is  expected  that  many  joints  (e.g.  between 
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wing  and  pylons,  between  riveted  panels,  between  bolted  pieces  of  spars,  ribs,  etc.)  would 
likely  participate  in  the  damping  effects.  However,  the  modeling  of  these  different 
contacts  would  not  only  be  a  tremendous  task  but  its  complexity  would  render  it 
extremely  difficult  to  determine  the  role  of  friction  in  LCO.  To  obtain  such  results,  it  is 
suggested  to  proceed  with  the  following  structural  simplifications  and  the  aerodynamic 
modeling  that  follows. 

(i)  The  friction  effects  will  be  assumed  to  originate  from  a  single,  rigid 
component  moving  relatively  with  respect  to  the  wing  and  to  be  characterized  by  a 
preload  N  or  applied  moment  M0  and  coefficients  of  friction  jus  and  . 

(ii)  The  changes  of  structural  geometry  that  occur  due  to  the  relative  sliding  at 
the  joints  will  be  neglected.  This  assumption  appears  justified  as  all  aircraft  components 
are  subjected  to  a  large  preload  so  that  a  large  amount  of  energy  can  be  dissipated  with 
only  small  relative  motions  at  the  joints.  Accordingly,  there  will  be  no  effect  of  the 
sliding  on  the  aerodynamics. 

Combining  the  above  two  assumptions  suggests  that  the  friction  mechanism 
should  appear  internal  to  the  wing,  e.g.  modeled  through  a  block  sliding  inside  the 
structure  in  a  rough  track  (Fig.  1.2)  or  through  a  rigid  component  experiencing  rubbing 
while  rotating,  e.g.  as  a  rotating  disk  squeezed  between  two  stationary  rings  (Fig.  1.3  and 
Fig.  1.4). 
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Both  of  these  models  will  be  used  in  the  ensuing  chapters  to  assess  the  role  of 
friction  on  aeroelastic  systems  exhibiting  either  a  pure  fluttering  behavior  (chapter  2)  or 
post  flutter  limit  cycle  oscillations  (chapter  3).  In  both  cases,  an  inviscid  incompressible 
aerodynamics  is  assumed  and  fully  coupled  aerodynamic/structural  dynamic 
computations  are  performed. 

It  should  be  noted  that  these  models  are  meant  to  be  generic,  i.e.  that  their 
equations  of  motion  may  be  considered  representative,  albeit  possibly  simplified,  of  the 
equations  of  motion  of  a  physical  wing/wing-store  system  expressed  in  modal 
coordinates  as  opposed  to  the  physical  coordinates  used  here.  A  one  to  one 
correspondence  between  the  properties  (mass  or  inertia,  stiffness,  ....)  of  the  friction 
device  and  of  a  physical  component  (store  or  missile  for  example)  is  not  suggested  here. 


CHAPTER  2 


FLAT  PLATE  WITH  INTERNAL  FRICTION 

As  a  first  step  forward  understanding  the  effect  of  friction  on  aeroelastic  system, 
the  simplest  airfoil,  i.e.  a  flat  plate,  will  be  considered  that  supports  a  torsional  friction 
device,  see  Fig.  1.3  and  1.4.  The  system  is  thus  characterized  by  the  three  degrees  of 
freedom:  plunging,  h,  (positive  downward  measured  at  the  elastic  axis),  pitching  about 
the  elastic  axis,  a ,  (positive  nose-up)  and  torsion  9  (positive  clock-wise). 


U 


Mean 

position 


Fig.  2.1  Flat  Plate  airfoil  having  three  degrees  of  freedom 


The  equations  of  motion  for  the  above  system  are 
Mh  +  Saa  +  khh  =  - L(t ) 


(2-1) 
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sah  +  Ia&  +  kca  +  kd(a-0)  =  Ma  (/)  +  Mf 

(2-2) 

Ide  +  kd{9-a)  =  -Mf 

(2-3) 

where  M  is  the  total  mass  of  the  system,  Sa  is  the  static  moment  of  the  flat  plate  about 
elastic  axis,  I a  is  the  moment  of  inertia  of  the  flat-plate  about  elastic  axis  and  I d  is  the 
moment  of  inertia  of  the  torsional  disk.  The  symbols  kh ,  ka  and  kd  denote  the  spring 
constants  of  the  plimging,  pitching,  and  torsional  motions,  respectively.  Further,  U  is  the 
free  stream  velocity  and  b  is  the  half  chord.  In  regards  to  forces,  L  and  Ma  denote  the  lift 

(positive  downward)  and  moment  (positive  nose-up)  respectively,  and  M,  is  the  friction 

moment  resulting  from  the  squeezing  of  the  rotating  disk  by  the  two  rings,  see  Fig  1.4. 
This  moment  can  be  expressed  as 


Mf  =  fJD  M0  sgn [O -a)  during  slip  (0*a) 

(2-4a) 

and  \Mf  |  <  /js  M0  during  stick  ( 9  =  a  ) 

(2-4b) 

where  /jdM0  is  the  moment  created  by  the  friction  induced  shear  stresses.  If  one  assumes 
a  uniform  stress  distribution  at  the  disk  ring  contact,  then  it  is  found  that 


M0=N-, 


R:-R 


R2-R? 


where  N  is  the  normal  force  between  the  disk  and  the  ring  and  Ri 


and  Re  are  the  internal  and  external  radii  of  the  rings.  Finally,  juD  and  jus  are  the 

dynamic  and  static  coefficients  of  friction,  respectively. 

Introducing  the  dimensionless  time  T-Uttb  in  Eq.  (2-1),  (2-2)  and  (2-3)  yields 


u 2  u 2 

M-fh'  +  S'—a'  +  W—m 

Sa  ^Th"  +  Ia  Ya"  +  k°a  +  ^  ~  ^  =  +  Mf 

Id~-e"  +  kd{d-a)  =  -Mf 

where  (•)  denotes  the  differentiation  with  respect  to  r  . 


(2-5) 

(2-6) 

(2-7) 


It  will  be  assumed  here  that  the  aerodynamics  can  be  assumed  as  inviscid  and 
incompressible.  The  resulting  aeroelastic  system  is  then  linear  in  the  absence  of  friction 
and  its  corresponding  response  can  be  sought  by  standard  Fourier  transfer  function 
concepts.  Specifically,  expressing  the  airfoil  motions  as 

h(  r)  =  h0e‘kr  (2-8a) 

and 


ikr 


a(r)  =  a0e 

where  the  reduced  frequency  k  is  defined  as 


k  _  cob 

~a 


(2-8b) 


(2-9) 


then,  it  can  be  shown  that  the  lift  and  moment  may  be  expressed  as 


L(r)e-*' _  ,2A 


pnbU 2  '  b 


=  -r  (-2-  -  ahaQ)  +  ika0  +  2  C(k) 


ao  +  T  kK  +  (t  ~  aA  )ika0 
0  L 


(2-10) 
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M{r)e,kt  ,1 


pjtb2U2  '2 


=  (-  +  ah)2C(k) 


aa+^-k\+{^-ah)ikaQ 
b  2 


~k\  ~  Wo)  -(~  ah)ika0  + y  a0 


(2-11) 


In  the  above  equations,  C(k )  denotes  the  Theodorsen  function9,10  simulating  the 
unsteadiness  of  the  wake  and  expressed  as 

H[2\k)  _  Kx{ik) 


C(k)  =  - 


(2-12) 


H\2\k)  +  iH(2\k)  Ka(ik)  +  K,Xik) 
where  //j2)(-)  and  (•)  denote  the  Hankel  function  of  the  second  kind  and  modified 


Bessel  function  both  of  order  j,  respectively. 

Adopting  a  rational  approximation  of  this  function  permits  the  reformulation  of 
Eq.  (2-10)  and  (2-11)  in  terms  of  ordinary  differential  equations  which  can  be  marched  in 
time  with  Eq.  (2-5),  (2-6)  and  (2-7)  to  obtain  the  response  of  the  coupled  aeroelastic 
system  both  without  (linear  analysis)  or  with  (nonlinear  analysis)  friction.  Specifically, 
the  following  approximation  of  C{k)  is  used  in  the  ensuing  analysis9,10 


C(k)=i-AiL-Ji± 


ik  +  5,  ik  +  s2 


(2-13) 


where  £  =  0.165,  s,  =0.0455,  £  =  0.335,  and  0.3. 

A  comparison  of  the  real  and  imaginary  parts  of  the  Theodorsen  function,  Eq.  (2- 
12),  and  its  rational  approximation,  Eq.  (2-13),  is  shown  in  Fig  2.2.  Clearly,  the  matching 
is  very  good  so  that  Eq.  (2-13)  is  appropriate  for  the  present  investigation. 
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Fig.  2.2  Comparison  of  the  exact  Theodorsen  function,  Eq.  (2-12),  with  its 
rational  approximation,  Eq.  (2-13) 


2.1  Time  Marching  Computation 

The  approximation  of  the  Theodorsen  function  by  a  rational  expression  is  only  the 
first  step  in  rewriting  Eq.  (2-10)  and  (2-1 1)  in  terms  of  h{r) ,  a(r)  and  their  derivatives. 
Indeed,  note  that  C(k)  does  not  appear  alone  in  Eq.  (2-10)  and  (2-1 1)  but  rather  through 
the  products  C(k)XQ  where  X0  denotes  either  h0  or  a0  and  that  there  variables  are, 


from  Eq.  (2-5),  the  Fourier  transform  of  h( r)  and  a(r) ,  i.e. 

K=ho(k)  =  F{h^)}  ;  «o  =a0(k)  =  F{cc(T)}  (2-14) 

Since  there  are  three  terms  on  the  right-hand-side  of  Eq.  (2-13),  one  can  write 
C(k)X0=F^'i  (r)  +  Tm  (r)  +  Ti+2  (r)} 


(2-15) 
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where  Tx(z)  arise  from  the  one  in  Eq.  (2-13)  while  TM (r)  and  Ti+2( z)  are  associated 
with  the  two  rational  terms  in  this  equation.  The  application  of  standard  rules  of  Fourier 


transformations  demonstrates  that 

T;  =  X(z) 

(2- 16a) 

Tm  :  tJx+sxTm=SxX' 

(2-1 6b) 

Ti+ 2  :  tJ+S2Tm=-Z2X' 

(2- 16c) 

There  are  thus  four  auxiliary  functions  of  time,  7j  (r)  and  T2  (r)  obtained  for  z= 0 
and  X(z)  =  h(z)  and  T3(z)  and  T4(z)  corresponding  to  i~2  and  X{z)  =  a(z)  . 

Other  noteworthy  Fourier  transform  properties  are 


ikC(k)X0  =F{t;+  T!4X  +  T!+2 }  (2-17*) 

ikX0=F{X'}  (2- 17b) 

and 

k  2X0=F{-X"}  (2- 17c) 

Combining  Eq.  (2-5),  (2-6),  (2-7),  and  (2-15)  finally  yields  the  equations  of  motion  in  the 
form 


M  1 
-  +  - 


pnV  b 


h"  + 


pnb 


3  ^ A 


a'  =  -a'  -  2  (a  +  T3  +  T4) -\(h'  +  T;+  7)') 

b 


-(1-2  ah)(a'  +  T;  +  W- 


pizbU 


(2-18) 
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pnb 


3  ah 


^  +  +  ^  +  \)a'  =  b(l  +  ^  )(<X  +  T%  +  T*  )  +  0  +  2°h  ){h'  +  T'+ 


A 


,1 


+b(--2al)(a'  +  T;  +  T;)-b(--ah)a' - '^a - + 

2  2  pnbU  pnbU 


M, 


pnbU 2 

(2-19) 


and 


U 


id-e'+kd{e-a)=-Mf 


(2-20) 


which  must  be  solved  together  with  Eq.  (2-16).  Note  during  stick  phases  that  6  and  a 
are  not  independent  variables  since  6-a.  Equation  (2-20)  then  permits  the  evaluation  of 
Mr. 


The  ensuing  numerical  study  was  conducted  at  various  speeds  U  and  for  different 


spring  stiffnesses  kd .  The  other  system  parameter  values  are  shown  in  Table  2.1. 
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parameter 

Explanation 

Value 

M 

Total  mass  of  the  airfoil  system 

20.0  kg/m 

Sa 

Static  moment  of  airfoil  about  the 

elastic  axis 

2  kg 

I 

Moment  of  inertia  of  airfoil  about 

the  elastic  axis 

1.1875  kgm 

L 

Moment  of  inertia  of  torsional  disk 

0.0625  kgm 

K 

Spring  constant  of  plunging 

direction 

1.8  N/m2 

K 

Spring  constant  of  pitching 

direction 

1.25N 

c r 

Semichord 

0.5m 

aH 

Nondimensional  distance  from  the 

midchord  to  the  elastic  axis 

-0.1 

Ms  (M) 

Static  coefficient  of  friction 

0.001 

Pd  (/*) 

Dynamic  coefficient  of  friction 

0.001 

P 

Density  of  air 

1. 27324% /m3 

M0 

Moment  term  in  Eq.  (2.4) 

1 

Table  2. 1  System  parameter  values  for  the  flat  plate 
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2.2  Linear  solutions  :  continuous  sticking  and  frictionless  slipping  behaviors 
Equation  (2-16)  and  (2-18)  through  (2-20)  represent  three  first  order  (two  during 
sticking  since  6  and  a  are  then  related  variables)  and  four  first  order  equation.  These 
equations  can  be  rewritten  as  a  set  of  10  first  order  equation  in  slipping  motions  and  8 
such  equations  in  sticking  phases.  These  equations  are  in  general  nonlinear  but  there  are 
two  linear  limiting  cases.  If  the  motions  are  such  that  sticking  occurs  continuously  (for 
example  in  the  limits  or  kd  — »  oo),  the  friction  nonlinearity  is  not  triggered  and 

the  equations  of  motion  are  linear.  This  sticking  solution  will  serve  as  the  baseline  in  the 
ensuing  analysis  as  it  corresponds  to  the  system  without  moving  parts. 

The  second  linear  solution  corresponds  to  the  continuous  slipping  motions 
occurring  with  /u  =  juD  =  jus  =  0 .  In  this  case  indeed,  the  nonlinear  moment  Mf  always 

vanishes. 

The  complexity  of  the  solutions  with  active  friction  warrants  the  detailed  analysis 
of  all  simple  cases  first.  In  this  regard,  note  that  the  equations  of  motion  for  both  of  the 
above  linear  problems  (continuous  sticking  and  frictionless  slipping)  can  be  rewritten  as 
y  +  Ay  =  0  (2-21) 

where  the  vector  y  has  the  10  components  h,  h ,  a ,  a ,  6 ,  0 ,  7],  T2,  and  7)  for  the 

frictionless  slipping  problem  and  only  the  8  components  h,  h ,  a ,  a ,  T2,  T3  and  7) 
for  the  continuous  sticking  solution.  Since  Eq.  (2-21)  is  linear,  no  LCO  is  possible  :  the 
response  y(t)  either  converges  to  0  or  diverges  independently  of  the  initial  conditions. 
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The  convergence  (stability)  or  divergence  (instability)  of  y{t)  is  determined  by  the 
eigenvalues^,  of  the  matrix^.  Specifically, 

the  system  is  stable,  i.e.  y{t)  converges  to  0  iff  the  real  parts  of  all  eigenvalues 
Xj  are>0. 

the  system  is  unstable,  i.e.  y{t)  diverges  iff  the  real  parts  of  any  eigenvalue  Xj  is 

<0. 

It  is  convenient  to  analyze  the  response  of  the  two  linear  cases  in  terms  of  natural 
frequencies  and  damping  ratios.  This  is  achieved  by  associating  to  each  eigenvalue  Xj  a 

natural  frequency  C0j  and  a  damping  ratio  .  through  the  linear  single-degree-of-freedom 
relation 

-Xj^-ojjCjtiOjfi^j  (2-22) 

or 

«HWI  <2-23> 

and 


(2-24) 


Stability  then  occurs  when  all  are  positive  and  instability  when  any  of  them  is 


negative. 
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This  stability  analysis  revealed  some  interesting  information,  most  notably  that 
the  system  in  ffictionless  slipping  mode  with  the  parameter  values  in  Table  2.1  is  less 
stable  than  its  continuous  sticking  counterpart  in  that  it  exhibits  its  first  negative  damping 
ratio  at  a  lower  flow  speed.  For  example,  the  first  negative  damping  ratio  arises  at 
approximately  £7=0.74  m/s  (continuous  slip)  and  £7=0.995  m/s  (continuous  stick)  for  kd  = 
0.0165  Nm/rad  but  they  occur  at  £7=0.99208  m/s  (continuous  slip)  and  £/=0.99232m/s 
(continuous  stick)  for  kd  =  0.985  Nm/rad.  It  must,  however,  be  recognized  that  there  are 
several  significant  differences  in  the  evolution  of  the  dimensionless  natural  frequencies 

(made  dimensionless  by  the  time  —;)  and  damping  ratios  corresponding  to  these  two 

values  of  kd .  For  the  larger  one,  see  Fig.  2.3,  classical  flutter  situations  are  encountered, 

two  modes  of  each  aeroelastic  system  (stuck  and  slipping)  appear  to  cross  as  the  smallest 
damping  ratio  becomes  negative.  Further,  note  that  the  natural  frequencies  and  damping 
ratios  of  the  sticking  and  slipping  systems  are  very  close  to  each  other  over  the  entire 
range  of  flow  speeds. 
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c) 


Zoomed  damping  ratios 

Fig.  2.3  Dimensionless  natural  frequencies  and  damping  ratios  of  the 
frictionless  slipping  (“si”)  and  continuous  sticking  (“st”)  systems  as  functions  of  the 
flow  speed,  kd  =0.985Nm/rad. 


For  the  smaller  value  of  kd ,  the  natural  frequency  of  the  torsional  device  alone  is 
the  middle  frequency  (as  opposed  to  the  largest  one  in  Fig.  2.3)  and  its  presence  affects 
much  more  significantly  the  response  of  the  frictionless  slipping  system.  This  system 
becomes  first  unstable  at  a  flow  speed  of  approximately  0.74m/s  but  this  does  not  occur 
with  a  mode  veering/crossing,  see  Fig.  2.4.  In  fact,  such  a  situation  does  not  take  place 
until  the  speed  reaches  0.995m/s.  Classical  flutter  instabilities  do  take  place  for  both  the 


2! 


a 


continuous  sticking  and  frictionless  slipping  systems  around  U  =  lm/s,  see  Fig.  2.5,  with 
the  sticking  system  instability  speed  lower  than  its  slipping  counterpart. 
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Fig.  2.4  Dimensionless  natural  frequencies  and  damping  ratios  as  functions 


of  the  flow  speed,  kd  =0.0165  Nm/rad. 


c) 


The  above  comments  demonstrate  that  there  is,  for  both  torsional  spring 
stiffnesses,  a  range  of  flow  speeds  for  which  the  continuously  sticking  system  is  stable 
while  the  one  in  frictionless  slip  mode  is  unstable.  Thus,  if  a  disturbance  is  applied  that 
does  not  induce  a  relative  motion  of  the  disk  with  respect  to  the  plate,  the  response  of  the 
system  will  decay.  On  the  contrary,  a  disturbance  that  includes  such  a  relative  motion 
may  increase  unless  the  system  gets  stuck  again.  These  observations  imply  that  the 
system  of  Fig.  2.1  might  exhibit  subcritical,  as  well  as  supercritical,  limit  cycle 


oscillations. 
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A  tentative  explanation  of  the  lower  flow  speed  at  which  instability  takes  place 
for  the  frictionless  slipping  system  as  compared  to  its  continuously  sticking  counterpart  is 
as  follows.  Note  first  that  the  stuck  configuration  (two  degree  of  freedom)  can  be  viewed 
as  the  limit  of  the  slipping  one  (three  degree  of  freedom)  in  which  the  connecting 
stiffness  of  the  device  kd  — >oc.  Accordingly,  the  two  natural  frequencies  of  the  stuck 

configuration  will  be  higher  than  the  two  lowest  natural  frequencies  of  the  slipping 
system.  Since  the  slipping  is  directly  coupled  to  the  pitching  dominated  motions  not  to 
their  plunging  counterparts,  it  is  expected  that  the  above  effect  would  be  more  significant 
on  the  pitching  frequency  than  on  the  plunging  one  and  the  former  would  thus  be 
decreased  by  allowing  relative  motion  of  the  disk. 

How  does  this  change  in  natural  frequency  affect  the  flutter  speed?  If  one  relies 
on  the  veering  of  the  natural  frequencies  at  flutter,  it  is  predicted  that  the  highest  natural 
frequency  of  the  coupled  aerodynamic-structural  system  decreases  while  the  lowest  one 
increases  as  the  flow  speed  is  increased  toward  the  flutter  point.  For  the  system  analyzed, 
the  pitching  frequency  is  the  highest  at  zero  flow  speed  and  thus  a  decrease  of  it,  through 
allowing  slipping,  may  result  in  a  lowering  of  the  flutter  speed.  This  explanation  is 
particularly  consistent  with  the  behavior  observed  for  kd  =  0.985  Nm/rad. 

The  above  discussion  is  the  characterization  of  the  system  behavior  in  linear 
frictionless  slipping  and  continuous  sticking  phases  and  provides  the  background  for  the 
analysis  of  the  effects  of  the  nonlinear  friction  force  to  be  accomplished  next. 
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2.3  Nonlinear  Responses 

An  important  aspect  of  the  numerical  integration  of  the  governing  equation  is  the 
accurate  capturing  of  the  times  at  which  a  zero  relative  speed,  6  -  a  =  0 ,  is  achieved.  It  is 
indeed  at  those  moments  that  the  force  of  friction  abruptly  changes  sign  (in  the  case  of 
continuously  slipping  motions)  but  also  that  slip  to  stick  transition  could  occur.  The 
appropriate  determination  of  the  stick  to  slip  transitions  is  also  important. 

An  interpolation  by  cubic  and  linear  polynomials  which  was  adopted  in  the 
ensuing  computations  is  presented  next. 

2.3.1  Transition  Handling 

If  a  transition  was  found  to  take  place  during  a  time  step  (i.e.  a  change  of  sign  of 
the  relative  velocity  or  the  friction  force  exceeding  the  sticking  threshold),  the  time 
/0  e  ,  tM  ]  at  which  it  occurred  was  estimated  by  interpolation  of  the  relative  velocity 

(using  a  cubic  polynomial)  or  the  friction  moment  (linear  interpolation). 

When  a  transition  (slip  to  slip,  slip  to  stick,  stick  to  slip)  occurs  during  the  time 
step  *e[/(.,t(+1],  the  estimate  of  the  response  cannot  be  obtained  directly  from  a  shooting 

of  the  solution  at  t-tl  because  the  friction  force  definition  and  the  equations  of  motion 

would  only  be  valid  before  the  transition  takes  place.  The  neglect  of  the  change  of 
character  during  the  time  step  would  produce  an  error  O(At)  and  thus  would  force  the 
selection  of  a  very  small  time  step  A t .  To  avoid  this  situation,  it  is  proposed  here  to  first 
march  the  solution  through  the  entire  time  step  as  if  no  transition  took  place,  then  to 
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estimate  the  transition  point  (  time  t0 )  within  the  time  step,  and  finally  to  recompute  the 
response  at  the  end  of  the  time  step  assuming  the  correct  state  in  the  interval  [t0  ,tM  ]. 

Consider  first  slip  to  slip  or  slip  to  stick  transitions.  The  transition  condition,  i.e. 
the  vanishing  of  the  relative  velocity,  is  based  on  the  response  variables  which  are 
continuously  differentiable.  Thus,  the  relative  displacement,  x(t),  can  be  approximated  by 
a  cubic  polynomial  of  time  in  the  interval  [  tt ,  tM  ],  i.e., 

x(t)  =  x0  +  aT  +  bl2  +  cV  (2-25) 

and 

x(t)  =  a  +  2bT  +  3cT2  (2-26) 

where  T  =  t  - tt  and  x0  =  x(t;) . 

The  constants  a,  b,  and  c  are  evaluated  from  the  values  of  the  relative 
displacements  and  velocities  at  the  beginning  and  end  of  the  time  step  assuming  that  no 
transition  takes  place.  Specifically,  when  t= 0  and  t=  At 

x(tM )  =  x,  =  x0  +  aAt  +  bAt2  +  cAt2  (2-27a) 

x(tt)  =  v0  =  a  (2-27b) 

x(tM  )  =  v,  =  a  +  2b  At  +  3c  At2  (2-27  c) 

Equation  (2-27b)  yields  the  value  of  a,  while  Eq.  (2-27a)  and  (2-27c)  lead  to  the 
values  of  the  remaining  coefficients  as 

(2-28a) 


bJ-rzl 

At 
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c-J-2/ 

At2 

(2-28b) 

where  /?  =  vt  -  v0  and 

*1  ~*0 

— - --n  =  r 

At 

(2-28c) 

On  the  basis  of  the  cubic  interpolation, 

it  is  possible  to  estimate  the  transition 

time.  Specifically,  from 

x(t0)  =  v0+2b70+3cT02  =0 

(2-29a) 

it  is  found  that 

-  -b±s]b2-3v0c 

3c 

(2-29b) 

Stick  to  slip  transitions  are  governed  by  the  value  of  the  force  of  friction  of  which 
only  the  value  is  known  (no  derivative).  It  will  thus  be  assumed  that  the  force  of  friction 
is  linear  during  the  interval  or 


1  A  t 


(2.30) 


where  Fx  and  F2  are  the  values  of  the  force  of  friction  at  the  beginning  and  end  of  the 
time  step  assuming  that  no  transition  took  place.  The  transition  occurs  when  F  reaches  its 
threshold  value  of  ±jusN  so  that 


■A/ 


(2.31) 
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Once  the  transition  time  has  been  estimated  from  Eq.  (2-29)  or  (2-31),  cubic 
interpolations  of  h(t),  a{t )  and  6{t)  are  carried  out  and  evaluated  at  the  appropriate  t0 
to  yield  approximate  values  of  the  response  at  the  transition. 

The  structural  computations  are  then  repeated  in  the  interval  /e[t0,t;+1]  with  the 
correct  state  (stick  or  slip). 

2.3.2  Time  Marching  Numerical  Result 

Based  on  the  analysis  of  Agelastos8,  it  was  of  particular  interest  to  determine  the 
type  of  response  (single  vs.  multiple  frequency,  continuous  slip  vs.  stick  slip  vs. 
continuous  stick)  and  the  largest  steady  state  responses,  both  as  functions  of  the  two 
fundamental  parameters:  the  flow  speed  U  which  dictates  the  level  of  instability  in  the 
system  and  the  spring  stiffness  kd  which  controls  the  transfer  of  energy  from  the  plate  to 

the  friction  device.  The  behavior  above  both  continuous  sticking  and  frictionless  slipping 
flutter  speeds  was  investigated  first,  and  Fig.  2.5  through  2.8  show  some  typical  examples 
of  observed  responses. 

The  simplest  motions  are  characterized  by  a  single  frequency  (SF)  and  continuous 
slip  (CS),  or  SFCS  for  short.  The  term  single  frequency  (see  Ref.  8)  refers  here  to  a 
solution  that  exhibits  a  single  fundamental  harmonic,  these  motions  are  thus  periodic  as 
can  be  seen  from  the  phase  plane  plot  of  Fig  2.5(d)  and  its  sampled  version  of  Fig.  2.5(e). 
Single  dots  in  this  latter  figure  indicate  that  the  crossing  of  the  zero  relative  velocity 
always  occurs  in  the  same  conditions.  Stick  slip  periodic  solutions  are  also  possible  and 
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are  denoted  as  SFSS  (single  frequency  stick  slip),  see  Fig.  2.6.  The  presence  of  the  stick 
phases  is  best  seen  in  the  plots  of  the  pitching  and  torsion  velocities,  see  Fig.  2.6(c);  they 
correspond  to  the  intervals  during  which  these  two  velocities  are  equal. 

Multiple  frequency  solutions  are  those  that  involve  at  least  two  frequencies  which 
are  not  rational  with  respect  to  each  other.  The  resulting  motions  are  thus  non-periodic 
and  exhibit  a  beating  as  seen  in  Fig  2.7  (continuous  slip)  and  2.8  (stick  slip).  As  might  be 
expected,  the  phase  plane  plots  of  these  multiple  frequency  (non-periodic)  motions  are 
much  more  complex  than  those  of  the  periodic  solutions.  In  fact,  the  MFSS  solution  of 
Fig  2.8  appears  from  Fig.  2.8(e)  to  be  aperiodic  which  its  MFCS  counterpart  would  seem 
to  be  chaotic,  see  Fig.  2.7(e). 

A  global  overview  of  the  response  of  the  system  can  be  obtained  by  analyzing  the 
two  most  essential  features  of  the  motions,  i.e.  the  type  of  solution  (MFSS,  MFCS,  SFCS, 
SFSS)  and  the  magnitude  of  the  response,  as  functions  of  the  two  critical  parameters  U 
and  kd .  The  concept  of  magnitude  is  not  clearly  defined  for  multiple  frequency  solutions, 

the  largest  noted  response  in  the  computed  steady  state  domain  will  be  used  instead.  This 
value  is  referred  to  as  the  largest  response  and  reduces  to  the  amplitude  of  response  for 
periodic  motions.  The  characteristics  of  the  response  and  its  largest  value  are  shown  vs. 
U  and  kd  in  Fig.  2.9(a)  and  2.9(b). 

There  figures  exhibit  many  of  the  features  presented  in  Ref.  8  Indeed,  it  is  seen 
that  the  largest  speed  at  which  a  stable  solution  exists  is  obtained  for  an  intermediate 
value  of  the  spring  stiffness  kd  (around  0.0165  Nm/rad),  that  a  transition  from  single  to 
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multiple  frequency  solutions  occurs  at  a  fixed  speed  as  the  spring  stiffness  kd  is 
increased,  and  that  both  continuous  slip  and  stick  slip  solutions  are  encountered. 

Considering  next  the  largest  pitching  response,  it  is  seen  in  particular  that  the 
response  level  decreases  first  monotonically  as  a  function  of  the  spring  stiffness  until  the 
transition  from  single  to  multiple  frequency  is  encountered.  For  larger  values  of  kd ,  the 
response  exhibits  beating  which  implies  the  very  sharp  increase  of  the  largest  response 
seen  in  Fig.  2.9(b)  near  0.0157  Nm/rad. 

There  are  a  few  differences  as  well  between  Fig.  2.9  and  their  counterparts  of  Ref. 
8.  In  particular.  Fig.  2.9(a)  and  2.9(b)  do  not  exhibit  the  symmetry  with  respect  to  the 
peak  stabilization  point  (kd&  0.0165  Nm/rad)  that  is  observed  in  Ref.  8. 

It  was  also  questioned  whether  the  “tuned  damper”  explanation  developed  in  Ref. 
8  in  connection  with  the  system  of  Fig.  1.1  could  also  be  used  to  justify  the  maximum 
stabilization  that  occurs  near  kd  «  0.0165  Nm/rad.  An  analysis  of  the  natural  frequencies 

of  the  system  in  slip  mode  revealed  that  they  were  indeed  close  to  each  other  but  with 
20%  difference  between  the  third  natural  frequency  and  the  two  that  are  veering/crossing. 
This  difference  is  too  large  for  the  tuned  damper  arguments  to  be  the  sole  justification  for 
the  optimality  of  kd  =  0.0165  Nm/rad  but  it  is  small  enough  for  these  arguments  to  be 
considered  as  the  basis  for  this  property. 

It  was  next  desired  to  investigate  the  behavior  of  the  response  in  situations  where 
the  stuck  system  is  stable  while  its  slipping  counterpart  is  unstable.  The  analysis  of  the 
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corresponding  response  types  and  largest  responses  was  achieved  first  around  kd  a  0.985 

Nm/rad.  It  was  found  that  limit  cycle  oscillations  do  occur  in  the  small  range  of  flow 
speeds  U  e  [0.922,  0.923]  m/s,  i.e.  from  the  flutter  speed  in  slip  mode  to  slightly  above 
the  flutter  speed  in  stick  mode.  All  finite  amplitude  limit  cycle  oscillations  found  were 
multiple  frequency  stick-slip  solutions  with  largest  pitching  response  of  the  order  of  0.2, 
i.e.  approximately  twice  as  large  as  those  seen  in  Fig.  2.9.  The  limited  zone  of  flow 
speeds  at  which  limit  cycle  oscillations  could  be  achieved  for  kd  a  0.985  Nm/rad  and  the 

larger  response  they  are  associated  with,  both  as  compared  to  kd~  0.0165  Nm/rad,  appear 
to  reinforce  the,  approximate,  validity  of  the  tuned  damper  assumption  as  the  third  natural 
frequency  in  slip  for  kd  ~  0.985  Nm/rad  is  far  from  the  other  two,  approximately  7  times 


larger. 


Time  history  of  plunging 


Time  history  of  pitching  and  torsion 


Time  history  of  velocity  in  pitching  and  torsion 
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d)  Phase  plane  plot  for  pitching  and  torsion 
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e)  Sampled  phase  plane  plot  at  the  moment  when  the  relative  velocity  becomes  zero 
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f)  Time  history  of  lift 


g)  Time  history  of  moment 


Fig.  2.6  n  -  0.001,  #<r=0.0155Nm/rad,  f/=1.0m/s  (SFSS) 
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phase  plane  plot  at  the  moment  when  the  relative  velocity  becomes  zero 
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f)  Time  history  of  lift 


g)  Time  history  of  moment 

Fig.  2.7  ju  =  0.001,  ftd=0.0155Nm/rad,  U= 1.0265m/s  (MFCS) 
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e)  Sampled  phase  plane  plot  at  the  moment  when  the  relative  velocity  becomes  zero 
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f)  Time  history  of  lift 


g)  Time  history  of  moment 

Fig.  2.8  n  =  0.001,  A,r=0.0162Nm/rad,  t/=1.013m/s  (MFSS) 
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a)  Characterization  of  the  stable  solutions 
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Fig.  2.9  Character  and  largest  responses  as  functions  of  U  and  kd 
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The  search  for  subcritical  limit  cycle  oscillation,  i.e.  occurring  at  speeds  lower 
than  the  flutter  speed  of  the  stuck  system,  was  continued  for  kd  =0.016  Nm/rad  in  the 
range  U=0.72  to  l.Olm/s.  In  fact,  it  was  found  that  an  unstable  LCO  does  exit  for 
U  <  Uj,  sluck  and  this  solution  is  a  separatrix,  i.e.  it  separates  the  motions  decaying  to  zero 

from  those  rapidly  growing  to  infinity.  In  general,  the  instability  of  this  separatrix  renders 
its  numerical  estimation  very  difficult.  It  was  nevertheless  estimated  here  by  varying  the 
magnitude  of  the  initial  condition  until  a  very  slow  decay  (at  first,  at  least)  was  obtained, 
the  system  response  during  that  phase  then  clearly  approximates  the  unstable  LCO.  On 
that  basis,  it  was  observed  that  the  motions  on  the  separatrix  would  be  characterized  as 
single  frequency  continuous  slip. 

For  speeds  larger  than  the  flutter  speed  of  the  stuck  system,  stable  single 
frequency  LCOs  were  obtained  for  U  in  the  range  [0.72,  1.013].  For  larger  values  of  U, 
multiple  frequency  solutions  were  observed,  see  Fig.  2.9.  Both  stable  and  unstable  LCOs 
pitching  amplitudes  are  shown  in  Fig.  2.10.  While  the  presence  of  the  separatrix  was  not 
confirmed  in  the  range  £/e[l,  1.013],  it  is  conjectured  that  it  is  present  and  that  the 
amplitude  of  these  unstable  motions  follows  the  dash  line  of  Fig.  2.10.  More  specifically, 
it  is  suggested  that  the  separatrix  meets  the  stable  single  frequency  LCO  branch  at  the 
flow  speed  where  this  solution  become  unstable  as  indicated  in  Fig.  2.10.  The  merging  of 
stable  and  unstable  branches  has  been  observed  in  a  1.5  degree  of  freedom  model.6,7 
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Fig.  2.10  Pitching  response  as  a  function  of  speed,  kd  =0.0160  Nm/rad 


2.3.3  Harmonic  Balance  Method 

The  results  presented  in  the  previous  section  have  been  obtained  by  marching  in 
time  the  system  response  from  a  set  of  initial  conditions  until  steady  state.  It  would  be 
desirable  to  have  a  procedure  that  allows  the  determination  of  the  steady  state  solution 
without  considering  the  transient  motions.  Two  such  approaches  are  possible. 

An  exact  approach  for  the  determination  of  the  steady  state  solution  can  be 
formulated  along  the  lines  suggested  by  Den  Hartog11  because  the  system  is  linear  except 
for  the  transition.  That  is,  a  set  of  yet  unknown  amplitudes  are  assumed  at  a  time  t= 0 
(conveniently  selected  as  the  time  at  which  the  relative  velocity  6 -a-  0)  and  the 
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response  is  analytically  computed  for  a  half-cycle  assuming  the  existence  of  zero,  one, 
two....  sticking  phases.  Then,  a  periodicity  condition  is  enforced  on  positions  and 
velocities.  For  a  continuous  slip  solution,  this  process  would  lead  to  10  equations  for  10 
unknowns  constituted  of  9  positions  and  velocities  at  t= 0  and  the  period.  While  these 
equations  are  linear  in  the  9  unknown  initial  conditions,  they  are  extremely  nonlinear 
with  respect  to  the  period  and  an  iterative  approach  is  necessary. 

A  simpler  but  approximate  method  to  obtain  the  steady  state  solution  is  to  rely  on 
the  harmonic  balance  method.  This  approach  rests  on  the  approximation  of  the  response 
as  a  Fourier  series  of  a  yet  unknown  fundamental  frequency.  Inserting  this  representation 
in  the  equations  of  motion  and  satisfying  them  to  the  highest  possible  harmonic  leads  to  a 
set  of  equations  for  the  unknown  fundamental  frequency  and  the  Fourier  coefficients. 
This  approach  is  most  reliable  when  the  response  is  closest  to  a  series  of  sine  and  cosine 
terms,  i.e.  in  continuous  slip  situations,  and  it  is  only  in  these  situations  that  it  will  be 
investigated  here. 

Considering  only  a  single  harmonic,  the  harmonic  balance  approximation  of  the 
response  is 

h{t)  =  A  cos  kz  +  B  sin  kz  =  Re^(/4-?5)e""  J  (2-32) 

a(t)  =  C  cos  kz  +  D  sin  kz  =  Re[{C-iD)eikT~\  (2-33) 

6{t)-a{t)  =  E  cos  kz  (2-34) 

Note  in  the  above  equation  that  6{t)-a(t)  does  not  contain  both  sin  Ar  and 
cos  kz  terms.  It  is  indeed  allowed  to  set  one  of  the  Fourier  coefficients  to  zero  in  free 
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response  problems  of  time  invariant  systems.  This  assumption  is  in  fact  equivalent  to 
imposing  t= 0  as  a  time  at  which  the  relative  velocity  vanishes. 

Then,  from  Eq.  (2-10)  and  (2-1 1)  the  lift  and  moment  can  be  expressed  as 
L  =  (Lhr  +  iLJ(A  -  iB)eik'  +  (Lar  +  iLai){C  -  iD)e'h  (2-35) 

M  =  (Mhr  +  iMhi)(A  -  iB)e~ikT  +  (Mar  +  iMai)(C  -  iD)eik*  (2-36) 

A  one  term  Fourier  series  approximation  of  the  moment  due  to  friction  is  obtained  as 


Mf  =sgn(8-a)juMQ  ~  --—sin  cat 


(2-37) 


n 


Substituting  the  above  expressions  for  all  displacements  and  friction  moment  into 
Eq.  (2.5)  through  (2.7)  leads  to  the  following  system  of  equations 
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(2-38) 


where  H  is  the  following  6x5  array 


H  = 
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(2-39) 


Proceeding  with  linear  combination  of  rows,  it  is  possible  to  rewrite  Eq.  (2-39)  as 


the  solution  pbtained  by  the  harmonic  balance  linearly  scales  with  the  coefficient  of 
friction.  In  fact,  this  result  is  also  valid  for  the  exact  solution  following  arguments 
developed  in  Ref.  8  and  it  is  sufficient  to  conduct  the  analysis  for  a  single  value  of  ju  as 
done  here. 
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b)  Maximum  pitching  amplitude 

Fig.  2.11  Comparison  of  LCO  frequency  and  pitching  amplitude  obtained  by 
time  marching  (TM)  and  the  harmonic  balance  method  (HBM),  kd  =0.0155Nm/rad 

The  above  procedure  was  applied  for  kd  =0.0155  Nm/rad  and 

U  e.  [1.0145,1.0265]  m/s.  Shown  in  Fig.  2.11  are  the  frequency  and  amplitudes  of 
response  obtained  by  both  the  harmonic  balance  method  and  the  time  marching  approach. 
Clearly,  the  matching  is  excellent  for  single  frequency  continuous  slip  solutions.  Figure 
2.12  provides  a  correspondence  between  the  actual  time  histories  for  =0.01 55  Nm/rad 

and  l/=1.0165m/s  obtained  by  the  time  marching  approach  and  computed  from  Eq.  (2-32) 
through  (2-34)  with  the  Fourier  coefficients  A,B,C,D,E  estimated  for  Eq.  (2-42)  and  (2- 
43).  This  figure  demonstrates  that  the  harmonic  balance  approximation  represents  well 
the  solution  as  a  function  of  time,  not  only  in  magnitudes. 
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c)  Torsion 


d)  Lift 
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Moment 


e) 


Fig.  2.12  Comparison  of  time  histories  obtained  by  the  time  marching 
approach  (TM)  and  the  harmonic  balance  method  (HBM). 
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CHAPTER 3 

AIRFOIL  WITH  INTERNAL  FRICTION 

The  previous  chapter  has  focused  on  the  assessment  of  friction  in  a  system  that 
does  not  exhibit  a  post-flutter  behavior  in  its  absence.  It  is  desired  next  to  evaluate  the 
potential  benefits  of  friction  in  aeroelastic  systems  that  do  exhibit  aero-driven  cycle 
oscillation.  The  work  of  Jadic  et  al.12  has  demonstrated  that  a  NACA0012  airfoil  in  an 
incompressible  inviscid  flow  may  exhibit  LCO  behavior  in  a  very  narrow  zone  of  speeds 
past  the  flutter  point.  This  example  will  thus  serve  here  to  assess  the  role  of  friction  in 
existing  LCO  and  both  the  sliding  and  torsional  friction  models  of  Fig.  1.2  and  1.3  will  be 
investigated.  The  aerodynamic  formulation  employed  here  is  identical  to  the  one  used  by 
Jadic  et  al.12,  Yao  and  Liu,13  and  Yao  et  al.14,15  and  is  briefly  reviewed  below  for 
completeness. 

3.1.  Aerodynamic  formulation 

The  flow  around  the  moving  airfoil  was  assumed  to  be  two-dimensional, 
unsteady,  inviscid  and  incompressible.  Accordingly,  the  velocity  field  can  be  expressed 
as 

£  =  +  (3-1) 

where  ^  is  a  pontential  satisfying  the  Laplace  equation 


A^  =  0 


(3-2) 
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In  addition  to  verifying  Eq.  (3-2),  the  potential  <j>  must  also  satisfy  the  following 
boundary  conditions 

(1)  The  flow  tangency  condition 

The  flow  is  assumed  inviscid  and  always  attached  on  the  surface  of  the  airfoil  so 

that 

(F0  -  Vs )  •  n  =  0  everywhere  on  the  airfoil  surface  (3-3) 

where  F,  is  the  local  velocity  of  the  airfoil  surface,  and  F0  is  the  local  flow  velocity  on 
the  airfoil,  i.e.  convection  velocity,  and  n  is  the  local  normal  to  the  airfoil  surface. 

(2)  The  Kelvin  Theorem 

The  total  circulation  T(t)  of  the  confined  flowfield  is  conserved 


(3)  The  far  field  boundary  condition 

VQ=U„  as  r->oo  (3-5) 

(4)  The  unsteady  Kutta  condition 

No  pressure  jump  must  take  place  across  the  wake  at  the  airfoil  trailing  edge,  or 
A Cp  =  0  at  trailing  edge  (3-6) 

The  current  method  adopts  the  linearized  scheme  developed  by  Kim  and  Mook16. 
The  determination  of  the  potential  <f>  and  velocity  field  V  <j>  satisfying  the  above 
equations  was  achieved  through  the  introduction  of  source  and  vorticity  distributions  on 
the  surface  of  the  panel  represented  by  high-order  elements  (cubic  spline  curved  panels 
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with  continuous  slopes  and  curvatures  at  the  joint  points,  see  Fig  3.1)  and  of  vortices  shed 
from  the  trailing  edge,  see  Fig.  3.2.  Once  shed,  the  vortices  were  assumed  to  convect  with 
the  local  flow  without  dissipating  thus  forming  an  unsteady  airfoil-free  wake  system. 

While  the  strengths  of  the  sources  were  specified  by  the  slope  of  the  airfoil,  those 
of  the  vortices  on  the  panel  (including  the  one  being  shed)  were  determined  at  each  time 
step  to  satisfy  Eq.  (3-3)  to  (3-6).  More  specifically,  the  flow  tangency  condition  (3-3)  led 
to  the  system  of  equations 

Ay  =  C  (3-7) 

where  A  is  the  aerodynamic  influence  coefficient  (AIC)  matrix  for  the  vorticity 
distributions  y  =  [yA  YTE  .  Further,  yA  is  the  vector  of  vorticity  strengths  on  the  airfoil 
section  and  Tte  is  the  vortex  shed  to  the  wake  at  time  t.  Finally,  C  denotes  the  vector  of 

net  induced  normal  velocities  at  the  control  points  due  to  the  sources  distributions,  wake, 
vortices,  and  airfoil  motions.  Finally,  the  integration  of  the  pressure  on  the  airfoil  results 
in  the  aerodynamic  lift  and  moments,  L  and  Ma  .  Further  details  concerning  this 
formulation  can  be  found  in  Yao  and  Liu13  and  Yao  et  al.14,15. 
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Cubic  Spline  Curved  Panel  with 
Continuous  Slopes  and  Curvatures 
at  End  Points 


Fig.  3.1  Curved  panel  and  linear  source  and  vorticity  distributions 


Source  Distribution  cr  Trailing  Edge 

Vorticity  Distribution  Y  Vortex  Fte 


Fig.  3-2  Aerodynamic  model 
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3.2  Airfoil  with  sliding  block 
3.2.1  Equations  of  motion 

The  structural  dynamics  model  used  in  the  present  numerical  experiment  assumes 
that  the  airfoil  exhibits  only  rigid  body  motion  in  pitching  (rotation)  and  plunging 
(vertical  movement).  In  the  absence  of  friction,  the  equations  of  motion  of  this  two- 
degree-of-freedom  system  are 

Mh  -Md  cos  aa  +  Md  sin  aa2  +khh  =  L  (3-8) 

Iaa -Md  cos  ah +  kaa  =  Ma  (3-9) 

where  h  is  the  plunging  displacement  of  the  airfoil,  defined  positive  upwards,  and  a  is 
the  pitching  displacement  of  the  airfoil,  defined  positive  nose  up.  The  coefficient  d  is  the 
distance  between  the  elastic  axis  and  the  center  of  mass,  M  is  the  mass  of  the  airfoil 
section  with  a  unit  spanwise  length,  Ia  is  the  corresponding  mass  moment  of  inertia  of 

the  airfoil  section.  Finally,  kh  and  ka  are  the  modal  stiffnesses  corresponding  to  the 

plunging  and  pitching  motions,  respectively. 

When  an  internal  sliding  block  is  present,  the  equations  of  motion  become 

Mh  -  Md  cos  ad  +  Md  sin  ad2  +  k  ,h -sm(  8 -a)(Tf  +kfr ) 

r  ••  2l  (3-10) 

-mcos(9-a)\-hcos(0-a)  +  rd  +  2rd-aa  I  =  L 

Md  cos  ah  -  /  a  -  k„a  -  a(Tr  +  k,r) 

r  ~  (3-11) 

-rm  -h cos (9 -a)  +  ra  +  2ra -ad  I  =  -Ma 


for  the  airfoil,  and 

m  sin(<9  -  a)h  -  mad  +  mr  -  mr  a1  +{Tf  +  kfr)  =  0  (3-12) 

for  the  sliding  block.  In  these  equations,  0  is  the  fixed  angle  between  the  chord  line  and 
the  track  in  which  the  block  slides  and  a  is  the  shortest  distance  between  the  elastic  axis 
and  the  track.  Further,  m  is  the  mass  of  the  block,  r(t)  is  its  displacement,  and  kf  is  the 

stiffness  of  the  spring  connecting  the  block  to  the  airfoil.  Finally,  Tf  is  the  friction  force 

between  the  block  and  the  track.  It  is  defined  as 

Tf  -  jlidN  ■  sgn(r)  when  slip  occurs  (3-1 3a) 

and 

| ^/ 1  <  when  sticking  occurs  (3-1 3b) 

where  N  is  the  normal  force  and  jus  and  jiD  are  the  static  and  dynamic  coefficient  of 
friction,  respectively. 

The  integration  of  the  equations  of  motion  broadly  follows  the  scheme  proposed 
by  Jadic  et  al.12  but  modified  to  account  for  the  presence  and  effects  of  the  sliding  block. 
That  is,  during  slip  phases,  Eq.  (3-10),  (3-1 1)  and  (3-12)  can  be  rewritten  as 


where 


F,=R,-^-R,_-^a 
c3  c3 


F2=-R2*^-R,-^-Kh 

C,  ct 


F2  =  -m  sin(<9  -  a) 


+  ma\  -1 - 2_ 

l  Mx  j 

l  ^  ) 

-Tf+  ma2r 


C,  =  M  +  mcos2  (0  -  a) 


C2  =  Md  cos  a  +  mr  cos(#  -  a) 


C3  =  Ia+  mr 


i?,  =  -Md  sin  aa2  +  sin(<9  -  a)(Tf  +  kfr) 
+  m  cos(0  -  a)(2ra  -  ad 2  )  +  L 

R2  =  a(Tf  +  kfr )  +  mr(2ra  -  aa2)-  Ma 


During  sticking  phases,  it  is  not  the  position  of  the  block  which  is  an  independent 
variable  but  rather  the  force  of  friction.  Specifically,  imposing  r  =  r  =  0  in  Eq.  (3-12),  it  is 
found  that 


(Tf  +  kfr)  =  -wsin(#  -  a)h  +  mda  +  mra 2 


(3-15) 


and  further  that 


where 


C  2 

c 2 

=C6~f~ 

F4  =  i?3  +—M  -—k  a 

**  ■>  /-r  a  a 

'-'6  ^6 

Fs  -  Ma  +—R,  ~—khh 

j  a  s-i  j  s~,  ft 

C4 

C4=M  +  m 

C5  =  Md  cosa  +  mr  cos(#  -  a)  +  ma  sin(#  -  a) 

C6=  Ia+mr2  +ma2 

R3-L- \md  sin  a  +  ma  cos(#  -a)-mr  sin(#  -  a)]  a 1 

3.2.2  Numerical  solution  of  the  equations  of  motion 

The  numerical  implementation  of  the  above  structural  and  aerodynamic  models 
necessitates  the  careful  handling  of  the  fluid-structure  interaction  and  of  the  friction 
related  transitions,  i.e.  slip  to  stick,  stick  to  slip  and  slip  to  slip. 

In  regard  to  the  first  issue,  a  full  aero-structure  feedback  was  desired  and  was 
accomplished  in  an  iterative  manner  within  each  time  step  of  the  computations12.  To 
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clarify  this  strategy,  let  ti  and  tM  denote  the  beginning  and  end  of  the  i'h  time  step.  Then, 
at  the  first  iteration  of  that  time  step,  an  estimate  of  the  structural  response  at  tM  was 
obtained  by  integrating  the  structural  dynamics  equations  from  t.  and  tM  assuming  right- 
hand-sides  of  Eq.  (3-14)  or  (3-16)  (including  the  aerodynamic  loading)  to  be  constant  for 
the  entire  time  step  and  with  values  obtained  from  time  t. .  This  estimate  of  the  structural 
response  was  then  used  to  update  the  boundary  conditions  for  the  flow  field  calculations 
which  in  turn  provided  an  estimate  of  the  aerodynamic  loading  at  tM .  At  the  second 
iteration,  these  values  and  their  counterparts  at  t,  were  combined  to  yield  a  linear 
interpolation  of  the  right-hand-sides  of  Eq.  (3-14)  or  (3-16)  from  which  the  structural 
response  at  /.+1  was  again  obtained.  The  second  iteration  was  completed  by  the  updating 

of  the  aerodynamic  loading  at  tM  and  a  third  iteration  was  then  started.  The  procedure 

was  continued  until  a  prescribed  precision  on  h,  a ,  and  r  (10-10  relative  error  between 
successive  iterations)  was  attained.  The  computation  at  the  next  time  step  was  then 
initiated  and  the  process  repeated  until  the  total  time  was  reached. 

The  handling  of  the  friction  induced  transitions  was  found  to  fit  very  naturally 
in  the  above  framework  and  was  achieved  by  the  cubic/linear  interpolation  strategy 


described  in  section  2.3.1. 


63 


3.2.3  Result  of  the  sliding  block  study 

The  analysis  of  the  system  of  Fig.  1.2  was  carried  out  for  m !  M  =  0.1.  Further, 
the  other  system  parameters  were  selected  so  that  the  system  in  sticking  mode  reduces  to 
the  one  considered  in  Ref.  12.  That  is,  m  +  M=  20.0  kg/m;  I  =  1.25  kg  m,  kh=  1.8  N/m2 

;  ka=  1.25N,  d  =  0.1  m,  c  =  1  m,  0  =  45 deg,  and  a=  0.  The  equations  of  motion  were 

nondimensionalized  by  using  the  chord  length  and  the  time  needed  for  the  flow  to  pass 
through  the  airfoil.  At  time  t=  0,  the  mass  m  was  assumed  to  be  at  the  elastic  axis  of  the 
airfoil.  Thus,  if  it  remained  there,  the  airfoil  properties  would  be  identical  to  those 
considered  by  Jadic  et  al.12  for  which  the  flutter  speed  was  found  to  be  Va  =  2.02 

(Fw  =-^-2-)  and  where  an  LCO  naturally  occurs.  The  system  was  then  analyzed  from 
coab 

Vn  =  2.02  to  2.4.  Additionally,  different  stiffnesses  of  the  internal  spring  were  selected 
that  yielded  natural  frequencies  of  the  secondary  system  close  to  that  of  the  LCO  in  the 
absence  of  friction,  i.e.  kf  =  0.36,  0.72,  1.00,  and  1.20  N/m.  For  all  cases  considered, 

small  initial  velocities  were  given  to  the  pitch  and  block  degrees  of  freedom  ( a  =0.1, 
r=  0.1  nondim).  Limit  cycle  oscillations  with  the  movable  friction  block  were  noted  for 
nondimensional  velocities  up  to  Va  =  2.2 . 

The  analysis  of  this  broad  set  of  data  demonstrated  the  primary  existence  of 
single  frequency  continuous  slip  and  continuous  stick  solutions,  see  Fig.  (3.3)  and  (3.4), 
with  the  stuck  solutions  occurring  at  the  highest  coefficients  of  friction  considered. 


pitching  Plunging 


Time  history  of  plunging 


Time  history  of  pitching 


Time  history  of  the  block  position 


c) 


d)  Time  history  of  the  relative  velocity  of  the  block 
Fig.  3-3  Response  of  the  system,  fx  =0.02,  ^/=1.2N/m  and  V=2.2  (SFCS) 


Time  history  of  plunging 


Time 


Time  history  of  pitching 
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Time  history  of  the  relative  velocity  of  the  block 

Fig.  3.4  Response  of  the  system,  fi  =0.08,  A/=1.2N/m  and  V=2.2  (Stuck) 

The  occurrence  of  LCO’s  in  stuck  configurations  of  the  block  was  not  initially 
expected  and  does  as  not  have  a  counterpart  in  Ref.  8.  This  surprising  result  is  due  to  the 
effect  of  the  block  location  on  the  inertia  of  the  entire  airfoil.  Indeed,  placing  the  block 
(in  a  stuck  configuration)  away  from  the  airfoil  elastic  axis  increases  the  inertia  of  the 
system,  lowers  its  pitching  natural  frequency,  and  increases  the  flutter  speed.  That  is, 
placing  the  mass  away  from  the  center,  without  any  motion  or  friction,  already  provides  a 
stabilization  of  the  system.  This  simple  observation  demonstrates  that  it  is  very  difficult 
to  assess  the  bonafide  effects  of  friction  on  LCO  amplitudes  for  the  system  of  Fig.  1.2. 
Accordingly,  no  further  discussions  of  this  model  will  take  place  and  the  airfoil  with  the 
torsional  friction  device  of  Fig.  1.3  and  1.4  will  be  considered  instead  as  it  does  not  suffer 
from  the  same  issues. 

3.3  Airfoil  with  rotating  disk 

Denoting  by  a  ( t ),  h(t),  and  6{t)  the  time  varying  pitching  and  plunging  of  the 
airfoil  and  the  rotation  of  the  internal  disk,  it  is  found  that  the  equations  of  motion  of  the 
system  are  (the  horizontal  translation  is  assumed  to  be  blocked) 

M  h-M  acosaa  +  M  asinaa2 +khh  =  L  (3.21) 

(/-  Id)cc  -Ma cos  a  h  +  (ka  +kd)a-kJ9  =  Ma  +Mj- 


(3.22) 


In  these  equations,  the  coefficient  a  is  the  distance  between  the  elastic  axis  and 
the  center  of  mass  (the  distance  PG ,  see  Fig.  1.3),  M  and  I  are  the  total  mass  and 
moment  of  inertia  of  the  airfoil  (including  the  disk),  Id  is  the  moment  of  inertia  of  the 
internal  disk,  and  kh,  ka ,  kd  are  the  stiffnesses  in  plunging,  in  pitching,  and  in  torsion  of 
the  disk.  Finally,  L  and  Ma  are  the  aerodynamic  lift  and  moment  acting  on  the  airfoil, 
and  Mf  is  the  moment  associated  with  friction,  see  Eq.  (2.4). 

3.3.1  Results  of  the  rotating  disk  study 

To  echo  the  study  of  the  lumped  mass  friction  system,  the  analysis  of  the  system 
of  Fig.  1.3  and  1.4  was  carried  out  for  Id  II  =  0.05 .  Further,  the  other  system  parameter 

were  set  to  m  =  20.0  kg/m;  I  +  Id=  1 .25  kg  m,  kh —  1 .8  N/m2;  ka  =  1 .25  N,  a  =  0. 1  m,  and 

c  =  1  m.  Thus,  the  behavior  of  the  airfoil  when  the  disk  is  stuck  is  identical  to  the  one 
obtained  by  Jadic  et  al.12  The  equations  of  motion  were  nondimensionalized  by  using  the 
chord  length  and  the  time  needed  for  the  flow  to  pass  through  the  airfoil.  The  system  was 
then  analyzed  right  at  the  flutter  speed  of  the  stuck  system  (nondim.  VM  =  2.02  ),  where  a 
LCO  naturally  occurs. 

From  the  time  history  of  the  response  in  stuck  mode,  it  was  determined  that  the 
disk  would  be  stuck  at  all  times  if  p  >  0.006 .  Accordingly,  the  analysis  of  the  airfoil  with 
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the  rotating  disk  was  conducted  for  jj.  <  0.006  The  stiffness  kd  was  varied  and  the  steady 
state  response  of  the  system  was  determined  by  numerically  integrating  the  equations  of 
motion  (structural  and  aerodynamic)  from  a  set  of  small  initial  conditions.  Shown  in  Fig. 
3. 3-3. 8  are  the  largest  responses  of  the  pitching,  plunging,  and  torsional  degrees-of- 
freedom  for  //  =  0.001  -  0.006.  Note  that  the  response  of  the  system  in  stuck  mode  (or 
without  the  disk)  is  also  shown  as  “h  stuck”  and  “alpha  stuck”,  see  Jadic  et  al.12 

It  is  seen  in  these  figures  that  the  responses  are  below  their  stuck  levels  for  low 
values  of  the  stiffness  kd  and  decrease  further  as  kd  is  increased.  A  minimum  of  both 

pitching  and  plunging  responses  is  obtained  at  an  intermediate  value  of  kd  but  it  is 
immediately  followed  by  a  sharp  rise  (especially  for  the  smaller  values  of  n )  of  all 
response  amplitudes  indicative  of  the  resonance  of  the  torsional  disk  system,  i.e.  when  its 
natural  frequency  is  very  close  to  the  flutter  frequency.  This  observation  is  in  complete 
agreement  with  the  “tuned  damper”  discussion  of  Agelastos8.  Increasing  further  the 
stiffness  kd  takes  the  system  out  of  resonance  and  all  three  responses  decrease  slowly 
toward  their  stuck  configuration  levels. 
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Figure  3.5.  Amplitude  of  response  in  plunging,  pitching,  and  torsion  (for  the 
disk)  as  a  function  of  the  stiffness  kd ,  =0.001 


Figure  3.6  Amplitude  of  response  in  plunging,  pitching,  and  torsion  (for  the 
disk)  as  a  function  of  the  stiffness  kd ,  n  =0.002 
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Fig.  3.7  Amplitude  of  response  in  plunging,  pitching,  and  torsion  (for  the 
disk)  as  a  function  of  the  stiffness  kd ,  // =0.003 
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kd 

Figure  3.8  Amplitude  of  response  in  plunging,  pitching,  and  torsion  (for  the 


disk)  as  a  function  of  the  stiffness  kd ,  n  =0.004 


Figure  3.9  Amplitude  of  response  in  plunging,  pitching,  and  torsion  (for  the 
disk)  as  a  function  of  the  stiffness  kd ,  fj.  =0.005 


Figure  3.10  Amplitude  of  response  in  plunging,  pitching,  and  torsion  (for  the 
disk)  as  a  function  of  the  stiffness  kd ,  /u  =0.006 
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To  condense  further  the  results  of  Fig.  3.5  through  3.10,  it  was  decided  to  record 
for  each  value  of  jj.  both  the  minimum  and  maximum  values  of  the  pitching  and 

plunging  amplitudes  over  the  domain  kd  e  [0,oo)  .  The  evolution  of  the  minimum 

amplitude  of  the  plunging  and  pitching  motions  with  increasing  coefficient  of  friction  is 
shown  in  Fig.  3.11.  It  is  seen  that  the  benefit  of  friction  is  most  significant  for  low 
coefficients  of  friction  as  the  corresponding  minimum  amplitudes  of  response  are  well 
below  the  stuck  values  from  Ref.  12.  But  the  penalty  for  this  choice  is  a  risk  of  very  high 
amplitudes,  see  Fig.  3.12  for  the  maximum  amplitudes,  as  the  minimum  and  maximum 
amplitudes  occur  within  a  very  narrow  range  of  kd  values. 

The  trend  of  the  minimum  amplitudes,  see  Fig.  3.11,  i.e.  a  monotonic  decrease 
with  decreasing  coefficient  of  friction  suggests  the  analysis  of  the  frictionless  system,  i.e. 
with  jjd=  /us=  0  .  Interestingly,  it  was  found  that  the  corresponding  airfoil  motions 

converge  to  large  amplitude  LCOs  when  kd  is  larger  than  0.01685  N/m  but  the  response 
appears  to  diverge  for  kd  <0.01685  N/m.  It  is  tentatively  suggested  that  this  behavior 

originates  from  a  difference  in  flutter  speeds  of  the  continuously  sticking  and  frictionless 
slipping  systems  (as  discussed  in  connection  with  the  flat  plate  in  Chapter  2)  and  the 
limited  domain  of  stable  LCO  of  the  NACA0012  airfoil  in  the  absence  of  friction. 

The  character  of  the  various  stable  limit  cycle  was  also  investigated,  see  Fig. 
3.13.  It  is  seen  that  single  frequency  stick  slip  solutions  e.g.  see  Fig.  3.14,  are  found  in  the 
most  of  the  domain  with  the  exception  of  the  values  of  kd  and  jj.  that  yield  large 
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torsional  responses.  In  that  zone,  the  motions  are  primarily  single  frequency  continuous 
slip  motion  e.g.  see  Fig.  3.15,  with  some  rare  multiple  frequency  solutions,  e.g.  see  Fig. 
3.16. 

This  characterization  is  expected:  when  a  low  torsional  response  is  induced, 
sticking  is  likely  to  occur  while  continuous  slip  motions  will  be  seen  for  large  rotations  of 
the  disk.  This  latter  situation  occurs  when  the  transfer  of  energy  from  the  airfoil  to  the 
disk  takes  place  efficiently,  i.e.  for  low  values  of  jj.  and  for  values  of  kd  that  induce  a 
torsional  natural  frequency  close  to  the  flutter  frequency.  Finally,  the  occurrence  of 
multiple  frequency  solution  at  the  center  of  the  single  frequency  continuous  slip  zone  is 
consistent  with  the  analysis  of  Agelatos8. 

The  pattern  of  the  vortices  forming  the  wake  was  finally  analyzed.  It  was  found 
that  these  patterns  clearly  parallel  the  steady  state  motions;  they  form  a  perfect  harmonic 
arrangement  for  single  frequency  solutions,  see  Fig.  3.15(c)  but  exhibit  beating  when  the 
system  response  is  dictated  by  multiple  frequencies,  see  Fig.  3.16(c). 


mu 


Time  histories  of  pitching  (a.dot)  and  torsional  (t.dot)  velocities 


c)  Position  of  wake  vortices  at  nondimensionlized  time  t=500 


Fig.  3.15:  SFCS  case,  /i  =0.001,  kd  =0.01125  Nm/rad 


Position  of  wake  vortices  at  nondimensionlized  time  f=500 


Fig.  3.16:  MFSS  case,  ju  =0.001,  kd  =0.0168  Nm/rad 
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CHAPTER  4 
CONCLUSIONS 

The  focus  of  this  thesis  has  been  on  the  assessment  of  the  effects  of  internal 
friction  on  the  response  of  aeroelastic  systems  exhibiting  either  explosive  flutter  or  limit 
cycle  oscillations. 

In  the  first  part  of  this  thesis,  a  flat-plate  airfoil  model  is  considered  which 
supports  a  torsional  friction  device  composed  of  a  disk  flexibly  connected  to  the  plate  by 
a  torsional  spring  and  squeezed  between  two  rough  surfaces.  The  behavior  of  this  system 
is  studied  when  placed  in  a  uniform,  inviscid  and  incompressible  flow.  A  stability 
analysis  of  the  system  with  the  disk  in  either  continuous  sticking  or  frictionless  slipping 
modes  was  carried  out  first  to  assess  the  expected  stability  domains.  It  was  shown  and 
justified  that  the  system  in  slip  mode  exhibits  instabilities  at  earlier  flow  velocities  than 
its  stuck  counterpart.  This  property  allows  for  the  existence  of  both  super-  and  subcritical 
limit  cycle  oscillations.  While  the  subcritical  limit  cycles  were  observed  to  be  unstable,  a 
zone  of  stable  supercritical  limit  cycle  oscillations  was  found  that  extends  about  3%  past 
the  flutter  speed  of  the  system  without  the  friction  device.  This  gain  shows  a  good 
stabilization  property  since  the  moment  of  inertia  of  the  selected  friction  device  system  is 
only  5%  of  the  moment  of  inertia  of  the  plate.  The  observed  limit  cycles  exhibit  either 
continuous  slip  or  stick  slip  behaviors  and  are  either  single  frequency  (periodic)  or 
multiple  frequency  (aperiodic  or  chaotic)  with  the  latter  ones  appearing  primarily  at  the 
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highest  flow  speeds  and  for  the  highest  frequencies  of  the  torsional  friction  device.  The 
above  results  were  obtained  by  time  marching  the  plate  equations  of  motions  with  a 
rational  approximation  of  the  Theodorsen  function  but  a  harmonic  balance  approach  was 
also  developed  that  led  to  very  good  approximations  of  the  single  frequency  continuous 
slip  limit  cycle  oscillations. 

The  second  part  of  this  thesis  focused  on  the  response  of  a  NACA0012  airfoil 
placed  in  a  uniform  inviscid  and  incompressible  air  flow  and  supporting  either  the  same 
frictional  device  as  the  flat  plate  or  a  block  sliding  in  a  rough  internal  track.  This  airfoil 
was  shown  in  earlier  studies  to  exhibit  limit  cycle  oscillations  at  and  slightly  above  its 
flutter  speed  without  any  frictional  device.  The  study  of  this  system  with  a  frictional 
device  provides  thus  a  basis  for  the  assessment  of  the  potential  effects  of  friction  on 
aeroelastic  systems  exhibiting  aero-driven  limit  cycles.  The  block  moving  in  an  internal 
track  was  first  considered  to  model  friction  in  the  system.  While  limit  cycle  oscillations 
were  observed,  it  was  also  shown  that  the  block  could  become  stuck  at  a  position  far  from 
its  original  one  and  thus  would  create  a  change  of  inertia  sufficiently  large  to  stabilize  the 
airfoil.  This  effect  does  not  involve  any  dissipation  due  to  friction  and  is  thus  not  relevant 
to  the  present  effort.  Accordingly,  this  frictional  model  was  not  considered  further  and 
was  replaced  by  its  torsional  counterpart  (as  in  the  flat  plate  analysis)  which  does  not 
suffer  from  the  same  defect.  The  results  of  time  marching  computations  demonstrate  that 
friction  can  substantially  decrease  the  level  of  the  limit  cycle  oscillations,  especially  with 
a  low  coefficient  of  friction,  but  that  increases  in  the  response  are  also  possible, 
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depending  on  the  selection  of  the  natural  frequency  of  the  torsional  friction  device.  As  in 
the  flat  plate,  continuous  slip  and  stick  slip  solutions  were  observed  most  of  which  were 
single  frequency  (periodic). 

The  results  of  this  study  demonstrate  that  friction  can  indeed  provide  a 
stabilization  of  an  impending  flutter  and  can  significantly  decrease  the  amplitude  of 
existing  limit  cycle  oscillations  of  aeroelastic  systems  with  an  appropriate  selection  of  the 
friction  device  parameters  most  notably  natural  frequency  and  coefficient  of  friction. 


85 


REFERENCES 


1  Denegri,  C.M.,  “Limit  Cycle  Oscillation  Flight  Test  Results  of  a  Fighter  with 
External  Stores,”  AIAA-2000- 1394 

2  Chen,  P.C.,  Sarhaddi,  D.  and  Liu,  D.D.,  “Limit  Cycle  Oscillations  Studies  of  a 
Fighter  with  External  Stores,”  AIAA-89-1727 

3  Edwards,  J.W.,  Schuster,  D.M.,  Spain,  C.V.,  Keller,  D.F.,  Moses,  R.W.,  “MAVRIV 
Flutter  Model  Transonic  Limit  Cycle  Oscillation  Test,”  AIAA-200 1-1291 

4  Cunningham  Jr.,  A.M.  and  Meijer.  J.J.,  “Semi-Empirical  Unsteady  Aerodynamics  for 
Modeling  Aircraft  Limit  Cycle  Oscillations  and  Other  Nonlinear  Aeroelastic  Problems,” 
International  Forum  on  Aeroelasticity  and  Structural  Dynamics,  Manchester,  U.K.,  Jun. 
1995. 

5  Mignolet,  M.P.,  Liu,  D.D.,  Chen,  P.C.,  “On  the  Nonlinear  Structural  Damping 
Mechanism  of  the  Wing/Store  Limit  Cycle  Oscillation,”  AIAA-99-1459 

6  Sinha,  A.,  and  Griffin,  J.H.,  “Friction  Damping  of  Flutter  in  Gas  Turbine  Engine 
Airfoils,”  Journal  of  Aircraft,  Vol.  20,  1983,  pp.  372-376. 

7  Sinha,  A.,  and  Griffin,  J.H.,  “Effects  of  Friction  Dampers  on  Aerodynamically 
Unstable  Rotor  Stages,”  AIAA  Journal,  Vol.  23,1 985,  pp.  262-270. 

8 

Agelastos,  A.M.,  Effects  of  Coulomb  Friction  on  Flutter  and  Limit  Cycle  Oscillations 
According  to  a  Structural  Dynamic  Model,  M.S.  Thesis,  Arizona  State  University,  May 
2005. 

9  Fung,  Y.C.,  “An  Introduction  to  the  Theory  of  Aeroelasticity”,  Dover  Publications, 

2002. 

10  Jones,  R.T.,  “The  Unsteady  Lift  of  a  Wing  of  Finite  Aspect  Ratio,”  NACA  Rept.  681, 
1940. 

11  Den  Hartog,  J.P.,  “Forced  Vibrations  With  Combined  Coulomb  and  Viscous 
Friction,”  Trans,  of  the  ASME,  APM-53-9,  1931,  pp.  107-115. 

12  Jadic,  I.,  So,  R.M.C.,  and  Mignolet,  M.P.,  “Analysis  of  Fluid-Structure  Interactions 
Using  a  Time  Marching  Technique,”  Journal  of  Fluids  and  Structures,  Vol.  12,  No.6,  pp. 
631-654,1998. 


13  Yao,  Z.X.,  and  Liu,  D.D.,  “Vortex  Dynamics  of  Blade-Blade  Interactions,”  AIAA 
Journal,  Vol.  36,  No.4  pp  497-504,  April,  1998 

14  Yao,  Z.X.,  Garcia-Fogeda,  P.,  Liu,  D.D.  and  Shen,  G.,  “Vortex/Wake  Flow  Studies 
for  Airfoils  in  Unsteady  Motions,”  AIAA  Paper  89-2225,  1989. 

15  Yao,  Z.X.,  Jadic  I„  So  R.M.C.,  Liu,  D.D.,  and  Tang  W.,  1995,  “Flow-Induced 
Vibrations  On  A  Turbine  Blade  In  a  Blade  Row  Using  A  Vortex  Dynamics  Model,” 
presented  at  the  international  forum  on  Aeroelasticity  and  Structural  Dynamics,  June  26- 
27,  Manchester,  United  Kingdom. 

16  Kim,  M.J.,  and  Mook,  D.T.,  “Application  of  Continuous  Vorticity  Panels  to  General 
Unsteady  Incompressible  Two-Dimensional  Lifting  Flows,”  Journal  of  Aircraft,  Vol.  23, 
pp.464-471,  1986 


